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1 Introduction 

This paper is primarily concerned with the relativistic generalizations of the iV-particle 
Calogero-Moser and Toda systems, both on the quantum and on the classical level. A 
survey of these systems can be found in [Rui94j . In addition, we briefly discuss the spe- 
cializations of our results to the nonrelativistic systems, surveyed in |OP81] and in |OP83] 
on the classical and quantum level, resp. 

The classical relativistic Calogero-Moser and Toda systems can be defined by Poisson 
commuting Hamiltonians of the form 

S h (x,p)= Yl V^HeMPPil k = l,...,N, (1.1) 

/c{i,...,JV} lei 

\i\=k 

with (3 = 1/mc, where m > is the particle rest mass and c > the speed of light. The 
elliptic version of the Calogero-Moser system describes N interacting particles on a line 
or ring, with Vj given by 

Vj{x) = Y[f(x m -x n ). (1.2) 

mel 
n<£I 

Here, the function / encoding the interaction is defined by 

m={ s[z+ :%;- p T ' (L3) 

where s(z) is essentially the Weierstrass sigma function (see (IA. 18[) in Appendix |A] for the 
precise relation). Throughout this paper we shall work with a fixed positive half-period 

u = 7c/{2r), r>0, (1.4) 
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and with two purely imaginary half-periods ia +J ia- in the relativistic quantum regime 
(cf. f l 1 . 2 3 f) below), whereas in the classical (h = 0) and nonrelativistic ((3 = 0) regimes we 
need only one purely imaginary half-period, parametrized as 

u' = ia/2, a>0. (1.5) 

With the 'coupling constant' p constrained by 

pei(0,a), (1.6) 

it follows that f(z) 2 is positive on the period interval z £ (0,n/r). By taking the posi- 
tive square root, we thus obtain well-defined positive coefficients V/(x) and Hamiltonians 
Sk(x,p) on the phase space 

Q = {{x,p) £ R 2N | x £ G}, (1.7) 

where G is the configuration space 

G = {x £ R N | x N < ■ ■ ■ < xi, x 1 - x N £ (0, 7r/r)}. (1.8) 

The classical relativistic version of the Calogero-Moser systems dates back to [RS86J. 
A quantization preserving commutativity was found in |Rui87j . For the elliptic systems 
it is given by the commuting analytic difference operators (henceforth AAOs) 

Sk{x) = 2j II f-( x m-x n ) Y[exp(-ihpd Xl ) f+(x m -x n ), k = 1, . . . , N, (1.9) 

ic{i,-,N}mei lei mei 

\I\=k n0 n0 

where h > is Planck's constant, and 

f ± (z) = (s(z±p)/s(z)) 1/2 . (1.10) 

With the above constraints on the parameters, they are formally self-adjoint. To promote 
them to commuting self-adjoint operators on the Hilbert space L 2 (G,dx), however, is an 
open problem. 

The unexpected existence of quite special kernel functions yields a novel perspective for 
solving this long-standing problem. This is explained in some detail in |Rui04] and [Rui09j . 
but in order to render our account more self-contained, we now digress to make the notion 
of 'kernel function' at issue in this paper more precise. 

Often, the term is loosely used for a function *ff(v,w) depending on variables v and 
w that may vary over spaces of different dimensions. Here, however, it is used in the 
restricted sense that the function at hand connects a pair of operators H\(v) and H2(w) 
acting on it via an equation of the form 

(H 1 (v)-H 2 (w))*(v,w) = 0, (1.11) 

which we refer to as a 'kernel identity'. Therefore, $?(v,w) can be viewed as a zero- 
eigenvalue eigenfunction of the difference operator Hi(v) — H 2 (w). In our setting, the 
operators Hi(v) and H 2 (w) will be either AAOs or PDOs. 

For the one- variable case v, w £ K, such kernel functions have a long history. For 
example, when Hi(v) and H2(w) are Mathieu operators (which can be viewed as Hamil- 
tonians for the reduced two- variable nonrelativistic periodic Toda system), a host of kernel 
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functions can be found in Section 4.1 of Arscott's monograph |Ars64] . Going yet further 
back, for certain pairs of Lame operators (which can be viewed as Hamiltonians for the re- 
duced two- variable nonrelativistic elliptic Calogero-Moser system), kernel functions occur 
in a paper by Whittaker |Whil5j . see also Section 23.6 of |WW35] . 

The appearance of kernel functions in a multi- variable setting is far more recent. 
For the nonrelativistic periodic Toda systems they can be found in a paper by Pasquier 
and Gaudin |PG92] , while for the nonrelativistic elliptic Calogero-Moser systems kernel 
functions emerged from Langmann's work on anyonic quantum field theory [LanOOj. We 
note that in the former paper explicit kernel identities of the above form do not yet appear 
(they are present implicitly), and in the latter only a pair of defining Hamiltonians was 
considered. 

Returning to the present paper, further work related to kernel functions will be cited 
in due course. Our starting point consists of the elliptic kernel functions from Ref. [Rui06j . 
They satisfy iV kernel identities of the form 

{S k {x) - S k (-y))V(x, y) = 0, k = 1, . . . , N, (1.12) 

and we shall recall their precise definition in Subsection 2.1. 

On the other hand, one key feature of the kernel functions should already be men- 
tioned now: Their building block is the elliptic gamma function (introduced and studied 
in |Rui97] ). and this function is symmetric under the interchange of the positive param- 
eters a and h/3 featuring in the A AOs Sk- This property is now often called 'modular 
invariance', and it entails that (II. 12ft also holds for the AAOs obtained by interchanging 
a and h/3. It is not hard to check that the latter AAOs commute with the previous ones, 
and so it is natural to insist on a joint diagonalization on L 2 (G, dx). 

In this paper, however, we are not directly concerned with joint eigenfunctions of the 
elliptic AAOs and their hyperbolic counterparts. Rather, we obtain some new insights 
concerning the hyperbolic kernel functions (as detailed in Subsection 2.2), and extend the 
theory of kernel functions to the relativistic Toda systems (in Subsections 2.3-2.5 and 
Appendix B). 

On the other hand, we do expect that the results of the present paper will be of piv- 
otal importance to solve open problems concerning the existence and properties of joint 
eigenfunctions of the AAOs of hyperbolic and Toda type. As a first example along these 
lines, we have shown in a recent paper how the hyperbolic kernel functions can be used 
for a recursive construction of modular invariant joint eigenfunctions of the commuting 
hyperbolic AAOs [HR12]. The idea that such a recursive construction might exist came 
from papers by Kharchev, Lebedev and Semenov-Tian-Shansky |KLS02j and by Gerasi- 
mov, Kharchev and Lebedev |GKL04j . which in turn were inspired by the pioneering work 
of Gutzwiller |Gut81j . As it turns out, the recursive constructions in |KLS02] and |GKL04] 
can be simplified and generalized in terms of kernel functions for the nonperiodic Toda 
systems and hyperbolic Calogero-Moser systems, all of which arise naturally within the 
context of the present paper. 

Accordingly, using the kernel functions detailed below as a starting point, we hope to 
arrive at a unified picture for the joint eigenfunctions of the various regimes. This will 
involve in particular detailed estimates on the kernel functions and their duals, which are 
necessary to control various analytical difficulties arising in the recursive constructions of 
the joint eigenfunctions and the study of their duality and Hilbert space properties. 



4 



In order to sketch our results on hyperbolic kernel functions, let us note first that the 
functions / (II. 3D and f± (11 . lOf) reduce to 



f(z) = (l-sinh 2 (7rp/a)/sinh 2 (7rz/a)) 1/2 , p6z(0,a), (1.13) 

f±(z) = (sinh(7r(2±p)/a)/sinh(7T2/a)) 1/2 . (1.14) 

It is also immediate that the elliptic kernel functions satisfying fll.l2j> have hyperbolic 
analogs: We need only replace the elliptic gamma function by its hyperbolic counterpart. 
The latter is modular invariant, like its elliptic generalization. (We review the relevant 
features of these gamma functions in Appendix A.) However, by contrast to the elliptic 
case, it is possible to obtain kernel functions relating the iV-particle AAOs Sk{x), k = 
1, . . . , N, to sums of (N— £)-particle AAOs, with £ — 1, . . . , N. This is because we can take 
UN, . . . , i/N-t+i to infinity, a procedure that has no elliptic analog. (For the trigonometric 
regime — which we do not consider — and for the case k = 1, a similar result was obtained 
first by Komori, Noumi and Shiraishi [KNS09J, among a host of other ones.) The kernel 
function for the case i = 1 is the key building block for the recursive construction of 
modular invariant joint eigenfunctions [HR12J. 

Turning to our Toda results, we first recall some general features. For the periodic 
and nonperiodic Toda systems the coefficients V/ in (II. ip are given by 



V l( X ) = EI fAXm+l-Xm) \ [ fT{Xm-Xm-l)- (1-15) 

The pair potential reads 



f T (z) = (l + 7 2 exp(27T2/a)) 1/2 , 7 G R, a > 0, (1.16) 
and the periodic and nonperiodic versions are encoded via the convention 

xo = xn, %n+i = %i, (periodic Toda), (1-17) 

xo = 00, xn+i = —00, (nonperiodic Toda). (1.18) 

In both cases the functions Sk(x,p) yield positive Poisson commuting Hamiltonians on 
the Toda phase space 

n T = {(x,p) 6l 2W }. (1.19) 

The classical relativistic Toda systems were introduced in |Rui90j . together with a 
quantization preserving commutativity. The latter is given by the commuting AAOs 



S k (x) = ^2 [ ] frixm+i - x m )Y[exp(-ih(3d Xl ) ] [ fa(x m - s m -i), 

JC{1,...,JV} m£l lei m£l 

\I\=k m+i£I m—\0 

k = l,...,N. (1.20) 

Since fx has period ia, the AAOs obtained by interchanging a and h/3 commute with 
the above ones, and so the question whether modular invariant eigenfunctions exist arises 
once again. We proceed with some remarks concerning this modular invariance feature, 
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which are also meant to elucidate the context of various g-Toda papers where modular 
symmetry is not at issue, cf. the paragraph below (11.241) . 

When the relativistic Toda systems were last surveyed in |Rui94] . modular invariance 
was not mentioned. Indeed, at that time no eigenfunctions were known at all, so it was not 
clear that one should be looking for eigenfunctions with this symmetry property. More- 
over, the Toda AAOs Sk given by (ll.20p are not even formally self-adjoint on L 2 (R Ar , dx) 
(as pointed out in Subsection 6.1 of |Rui94] ) . so that their Hilbert space status seemed 
quite opaque. 

It transpired from the above-mentioned work by Kharchev, Lebedev and Semenov- 
Tian-Shansky |KLS02] in the context of quantum group representation theory that a 
slight modification of the quantum coupling dependence remedies the lack of formal self- 
adjointness of the Toda AAOs. Moreover, the modular invariance property showed up in 
the eigenfunctions presented in [KLS02J and was tied in with Faddeev's notion of modular 
double of a quantum group |Fad99] . Likewise, van de Bult has shown that the modu- 
lar invariance of the 'relativistic' hypergeometric function (introduced in Subsection 6.3 
of [Rui94j ) can be understood from this quantum group perspective [vdB06j. 

Both formal self-adjointness and modular symmetry emerge naturally from our new 
results on Toda kernel functions. We expect that these results will be crucial to solve 
open problems concerning the modular invariant joint eigenfunctions obtained in [KLS02J, 
including duality properties, orthogonality and completeness. 

The change in coupling dependence entailing formal self-adjointness consists in replac- 
ing the classical pair potential in (jl.20p by the quantum counterpart 

1 /2 

f T {z)= (l + exp (^[20 + 277 + ^/3]] J . (1.21) 

This amounts to the replacement 

7 2 -> exp(7r(2?7 + ihf3)/a) (1.22) 

in (I1.16p . Clearly, this substitution does not change the AAO-commutativity features 
mentioned above. The parameter r\ G R plays the role of coupling constant, with the limit 
Tj — y — oo yielding the free theory. Moreover, in the classical limit H — > the shift into 
the complex plane disappears. Note in this connection that there is no classical analog 
of modular invariance. (More precisely, the Hamiltonians obtained from the functions 
Sk(x,p) by interchanging a and (3 do not Poisson commute with Si(x,p), . . . , Sn-i(x,p) 
for 1 < k < N and a ^ (3.) 

The new kernel functions ty(x, y) for the relativistic Toda systems obtained in Subsec- 
tions 2.3 and 2.4 for the periodic and nonperiodic case, resp., have the hyperbolic gamma 
function as their building block. Just as in the quantum elliptic and hyperbolic cases, we 
switch in the quantum Toda case to positive parameters 

a + = a, a_ = h/3, (1-23) 

in terms of which modular- invariant formulas involving the hyperbolic and elliptic gamma 
functions are more readily expressed. 

We obtain the Toda kernel functions in a somewhat tortuous way. Basically, we exploit 
the previous elliptic and hyperbolic Calogero-Moser results to arrive at them. The con- 
nection between the hyperbolic Calogero-Moser Hamiltonians and their nonperiodic Toda 
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counterparts was already observed and used by the second-named author in 1985, which 
led him to the relativistic Toda systems via the relativistic hyperbolic Calogero-Moser 
systems |Rui 90j . The relation between the defining Hamiltonian of the nonrelativistic 
elliptic Calogero-Moser system and its periodic Toda counterpart was first pointed out by 
Inozemtsev |Ino89] . and then generalized to all of the commuting relativistic Hamiltoni- 
ans, cf. the survey |Rui94j. 

Here we need these results as well, but to control the pertinent limit for the kernel 
functions seems not feasible on the elliptic level. Instead, we determine the limits of the 
functional equations expressing the kernel function property to obtain Toda functional 
equations. These can then be viewed as corresponding to kernel functions for periodic 
Toda AAOs related to the above ones by a similarity transformation. The details can be 
found in Subsection 2.3. 

In Subsection 2.4 we first show how nonperiodic Toda kernel functions can be obtained 
as a limit of the periodic ones. Just as in Subsection 2.2, we can also obtain kernel 
functions connecting iV-particle to M-particle AAOs, but here this seems only feasible 
for the case \N — M\ < 1. At the end of this subsection we detail how the previous 
nonperiodic results and a few new ones follow directly from their hyperbolic Calogero- 
Moser counterparts. Some readers might prefer this avenue, since it does not involve the 
periodic Toda and elliptic regimes. 

At this point we would like to mention that the limit transitions from Calogero-Moser 
to Toda type systems have recently become important from the viewpoint of quantum 
groups and Cherednik algebras, cf. the lecture notes [CM09] and various references given 
there. By contrast to our perspective and that of the paper [KLS02] cited earlier, this 
work involves a single deformation parameter q not on the unit circle. Here we are dealing 
with two parameters 

q + = exp(i7ra + /a_), g_ = exp(27ra_/a + ), (1-24) 

and for Hilbert space/quantum mechanical purposes it is of pivotal importance that a + 
and a_ be positive, so that \q±\ = 1. Moreover, the kernel functions and eigenfunctions 
are invariant under the interchange of a + and a_, which arises from the hyperbolic gamma 
function G(a + ,a_; z) featuring building block. 

For the case where the building block is the trigonometric gamma function (better 
known as the g-gamma function), the Hilbert space status of the g-Toda AAOs and their 
joint eigenfunctions is opaque, but in this case there are intimate connections to various 
issues in representation theory and algebraic geometry. Some early references include 
|Eti99j . |Sev00j . |OR02] . |GL03] . In particular, in a series of papers by Olshanetsky and 
Rogov (which can be traced from [OR02j ). the shift in the imaginary direction featuring 
in (11.201) was for the first time traded for a shift by h in the real direction, hence yielding 
a parameter q = exp(-h) < 1. Their work concerns the rank-1 (2-particle) case, whereas 
Etingof 's paper |Eti99] appears to be the first where the arbitrary-rank case is dealt with. 

The notion of 'dual relativistic Toda systems' at issue in Subsection 2.5 is not widely 
known. On the classical level these systems emerged from the explicit construction of an 
action-angle map for the nonperiodic Toda systems |Rui90j . They are integrable systems 
for which the actions p and angles x play the role of the positions and momenta in the 
original system, respectively. Specifically, the Poisson commuting Hamiltonians can be 
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chosen as 



ft(M> = ^S i^fej^l g "»(**■/«), * = 1,...,AT. (1.25) 



|/|=fc 



The classical hyperbolic relativistic Calogero-Moser systems are self-dual, since the 
action-angle map is essentially an involution [Rui88]. More specifically, the symmetric 
functions of the 'dual Lax matrix' 



A(x) = diag(exp(27rxi/a), . . . , exp(27TXAr/o;)), (1-26) 
turn into the dual Hamiltonians 
S k (p,x) = II (l-sinh 2 (7rp/«)/sinh 2 (/3(p m -p n )/2)) 1/2 J]exp(2vrx,/«). (1.27) 

IC{1,...,N} m€l lei 
\I\=k n 

The limit transition from the hyperbolic Calogero-Moser to the nonperiodic Toda systems 
has a counterpart for the classical and quantum duals. Just as for the original system, 
one needs to quantize the Hamiltonians (I1.27P via the analog of the reordering in ( II. 9p 
(cf. f ll.13p - p.14p ) before substituting x\ —> —ihdp r This yields again commuting hyper- 
bolic AAOs, and then the desired commuting AAOs for the dual Toda case follow from 
the pertinent limit. Thus we obtain quantum versions of (I1.25P and corresponding 
kernel functions that are once again built from the hyperbolic gamma function. 

Without a change in notation, however, the procedure just sketched would lead to 
awkward formulas. Indeed, when we switch to the parameters a± via (11.231) . then the 
factor sinh(/3(jS n — p m )/2) (for example) becomes sinh(a_(p n — p m )/2fr). Physically speak- 
ing, the unpleasant occurrence of h can be understood from the parameters a + and a_ 
having the dimension [position], whereas p has the dimension [momentum]. To obtain 
the desired dual modular symmetry, we should trade p for a dual variable afip/2'K with 
dimension [position]. We denote this new position by v , so that we need the substitution 

p = 2th;/ (a/3). (1.28) 

Using this variable, the dual hyperbolic AAOs again take the form (II. 9p with f± given 
by (H.14p and x replaced by v. 

With this change of notation in place, the dual Toda kernel functions are symmetric 
under the interchange of a+ and a_, so that they are also kernel functions for the modular 
transforms of the H k , obtained by interchanging a+ and a_. Somewhat surprisingly, for 
the dual nonperiodic Toda case we easily obtain kernel functions connecting the dual N- 
particle AAOs to their M-particle versions for any M < N, whereas we can only handle 
the M = N — 1 case in Subsection 2.4. 

Our results for the dual Toda case are collected in Subsection 2.5. A close relative 
of the kernel function connecting the iV-particle and (N — l)-particle dual AAOs has 
appeared in the above-mentioned work by Kharchev et al. [KLS02J. It is used in a re- 
cursive construction of joint eigenfunctions for the nonperiodic Toda AAOs, without a 
consideration of duality and Hilbert space properties. 

Before sketching the results of Section 3, we add an important remark concerning the 
AAOs and kernel functions at issue in Section 2. The AAOs are invariant when all of 



S 



their coordinates x n are shifted to x n + £. This entails that the kernel function property is 
preserved under such coordinate shifts. We make use of this freedom to choose convenient 
kernel functions. Another common feature is invariance of the kernel function property 
under multiplication by any function of the form 

N \ 

yj(^n — Dn) I , meromorphic. (1-29) 

vn=l / 

Hence, once we have identified one kernel function we immediately obtain an infinite- 
dimensional family of kernel functions. 

Section 3 is concerned with so-called Backhand transformations for the classical rela- 
tivistic Calogero-Moser and Toda systems. These are canonical transformations (x, p) i— > 
(y, q) that preserve the Poisson commuting Hamiltonians, derived from a generating func- 
tion F(x, y) via 

OF dF 

<■> nr/ ' ' V < 13 °> 

For the nonrelativistic Calogero-Moser systems such transformations appear to date back 
to work by Wojchiechowski Woj82|. For the nonrelativistic infinite Toda chain a Backhand 
transformation can already be found in Toda's monograph |Tod81j . It seems Gaudin was 
the first to realize that it can also be applied to the finite Toda systems, and that it can be 
tied in with the classical limit of a kernel function for their quantum versions, cf. Ch. 14 
in his monograph |Gau83] . Pasquier and Gaudin [PG92j then used the kernel function to 
study eigenvalues and eigenfunctions. 

For the nonrelativistic rational Calogero-Moser system a Backhand transformation 
was obtained via special solutions of the KP equation by Nijhoff and Pang [NP94J, 
|NP96] . They reinterpreted the generating function as a Lagrangian for a discrete map, 
which they viewed as a time-discretization of the defining Hamiltonian. In the same 
spirit, in Nijhoff/Ragnisco/Kuznetsov [NRK96] Backhand transformations (alias 'time- 
discretizations') for the relativistic Calogero-Moser systems were introduced and studied. 

Later on, Kuznetsov and Sklyanin elaborated on the general theory of Backhand trans- 
formations |KS98j . In particular, they reconsidered the relation between kernel functions 
on the quantum level and generating functions on the classical level. Gaudin already 
pointed out this relation in the special case of the nonrelativistic periodic Toda sys- 
tem |Gau83j , but in [K S98] it was suggested more generally that the semi-classical behav- 
ior of a kernel function $?(x, y) of the type we consider should be given by a formula of 
the form 

#(/*; x,y) ~ exp(-iF(x,y)/H), h^O, (1.31) 

where F(x,y) generates a Backlund transformation for the classical version via (jl.30p . 

Now it seems quite unlikely that this is generally true, as kernel functions exist in 
profusion. Indeed, assuming one has found an orthonormal base {0 n (x)}^L o of joint 
eigenfunctions for the commuting elliptic Hamiltonians (say), a function of the form 

oo 

K((a , a 1 ,...);x,y) = ^2 a n (f) n (x)(j) n (y), (1.32) 

n=0 

is a Hilbert-Schmidt kernel function for any (a ,ai, . . .) G £ 2 (N). Since the numbers a n 
are arbitrary, it is not even clear what one would mean by the semi-classical behavior of 
such a general kernel function. 
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We are, however, dealing with very special kernel functions, which can be expressed 
in terms of the elliptic gamma function and its specializations. In particular, there is a 
notion of 'classical limit' of the hyperbolic gamma function, which is tied to its appearance 
in the quantum scattering of the relativistic hyperbolic Calogero-Moser system. Indeed, 
in order to obtain the classical scattering (position shift) for h — » via a coherent state 
correspondence (which goes back to Hepp's fundamental paper |Hep74| ), a quite special 
limit is required. This is detailed in Eq. (4.73) of |Rui97] . and in terms of the hyperbolic 
gamma function it amounts to a certain zero step size limit, cf. Prop. III. 7 in |Rui97] . 
The elliptic counterparts of these limits are Eq. (4.98) and Prop. III. 13 in |Rui97] . (In 
Appendix A we have recalled the two pertinent limits, cf. (IA.29j) and (IA.14j) .) 

The point is now that with the associated ^-dependence in force, the asymptotic be- 
havior encoded in (jl.3ip does yield the generating function of a Backhand transformation, 
as we shall show in Section 3 for each of the different cases at issue. More is true: For the 
relativistic Calogero-Moser case these generating functions are basically the ones arrived 
at in |NRK96] . Also, the Backhand transformations for the relativistic Toda regimes can 
be tied in with results by Suris |Sur96] . and their nonrelativistic limits yield the ones 
already known from the papers cited earlier. 

There is however an unsettling phenomenon associated with these Backlund transfor- 
mations, which seems not to have been pointed out before: They correspond to Calogero- 
Moser and Toda systems of an unphysical nature, inasmuch as there seems to be no choice 
of parameters that yields complete flows and phase space coordinates that stay real for 
all times. For the Toda regimes this disease can be remedied by an analytic continua- 
tion, a state of affairs that was already noted and used by Gaudin in the nonrelativistic 
case [Gau83j . (More precisely, he starts from the Backlund transformation with the phys- 
ical positive coupling, and then finds an associated quantum kernel with the 'wrong' 
coupling; this can then be remedied by analytic continuation of positions. Since we start 
with a positive coupling on the quantum level, we need to reverse this procedure.) 

For the Calogero-Moser case, however, this is no option. Indeed, even for the very 
simplest degeneration, namely, the nonrelativistic rational N = 2 Calogero-Moser system, 
it seems impossible to avoid the 'negative coupling' behavior. For the discrete map at 
issue it shows up in real initial positions becoming complex after a number of discrete time 
steps. Even so, the circumstance that the pertinent quantum kernel functions give rise 
to generating functions of Backlund transformations is highly remarkable and deserves a 
further scrutiny from the viewpoint of global analysis. 

The relevant limits and their Backlund features yield somewhat unwieldy formulas, 
which is why we shall not detail them here. In Section 3 we reconsider successively the 
same regimes as in Section 2, omitting details whenever there is considerable similarity 
to previous cases. 

Section 4 is concerned with the nonrelativistic version of our results. Here, too, we 
leave out details when they can be readily supplied by specialization. Subsection 4.1 
deals with the nonrelativistic counterparts of the kernel functions of Section 2, whereas 
Subsection 4.2 is concerned with the nonrelativistic limits of the Backlund transformations 
of Section 3. As already mentioned, the kernel functions and Backlund transformations 
we arrive at in Section 4 are not new. On the other hand, we arrive at the relevant features 
in a novel way. 

In Appendix A we review some properties of the hyperbolic and elliptic gamma func- 
tions from |Rui97] we have occasion to use. In Appendix B we prove that the center- 
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of-mass relativistic periodic Toda kernel functions that connect the iV-particle AAOs 
give rise to Hilbert-Schmidt operators. We anticipate that this result will be the key to 
proving that the commuting AAOs can be promoted to commuting self-adjoint operators 
on L 2 (M. N ), a problem that has been wide open for several decades. 

To conclude this Introduction, a summary of our notation changes is in order. In- 
deed, as transpires from the above, we use several notational conventions that depend 
on the regime at issue. Our choices not only simplify formulas, but also reflect physical 
distinctions. Specifically, in all regimes we have a length scale a coming from the inter- 
action in the defining Hamiltonian. In the quantum relativistic cases, however, we have 
an additional length scale, namely h/3 (physically speaking, the Compton wave length of 
the particles under consideration). As explained above, modular symmetry interchanges 
these two parameters, which is why it is convenient to work with two equivalent length 
scales a± in Section 2, cf. (I1.23p . Now in the hyperbolic and nonperiodic Toda regimes we 
also have a notion of dual system, with the 'spectral variables' pi, . . . ,p N denoting asymp- 
totic momenta. The self-duality of the relativistic hyperbolic regime, however, makes it 
more natural to work with the dual position v defined by (11.281) . (This has in partic- 
ular the consequence that asymptotic plane waves do not have the usual dimensionless 
combination x ■ p/h in the exponent, but x ■ v/a + a- instead.) 

By contrast to Section 2, we study in Section 3 and Section 4 the classical {h = 0) 
and nonrelativistic (/3 = 0) settings, so that the length scale hj3 disappears. In the 
elliptic regime we therefore revert to the parameter a, whereas in the classical and in the 
nonrelativistic hyperbolic and Toda cases we trade a + = a for a parameter 

H = 2n/a, (1.33) 

with dimension [position] -1 . This change not only avoids a plethora of factors tt, but is 
also in accord with the self-duality of the classical relativistic hyperbolic regime. Indeed, \i 
is the parameter naturally dual to 0, as can already be gleaned by comparing the defining 
Hamiltonians (11. ip and the dual ones (ll.27p . (Cf. also the Lax matrix (I4.137P and dual 
Lax matrix (I4.139P to appreciate this self-duality feature.) 

2 Kernel functions 
2.1 The elliptic case 

In this subsection we review various elliptic quantities that play a role in our study of 
the periodic Toda case. As explained above, it is convenient to use notation that encodes 
modular invariance, cf. (I1.23p . To start with, we switch from the N commuting AAOs Sk 
given by (II. 9ft to the 2N commuting Hamiltonians 

H k ,s(x) = 22 EI fo-( x m ~ Xn ) II ex P(-ia-sdx m ) fs,+( x m ~ x n ), (2.1) 

IC{l,...,N} m£l mel m£l 

\I\=k nil nil 

where k = 1, . . . , N, 5 = +, —, and 

A±W = ( SS{Z ^ P) ) V2 , s s (z) = s(r, a s ; z). (2.2) 
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(See Appendix A for the definition and properties of the functions s±(z).) Next, we 
introduce 2N additional AAOs by setting 

H- k>s {x) = H k>s {-x), k = l,...,N, 6 = +,-. (2.3) 

Thus we have 



H- k ,s(x) = ^2 ] I fs,+{xm ~ x n ) \\ exp(ia^ s d Xm ) Y\ fs-(x m - x n ), (2.4) 
i 

and in particular 



IC{l,-,N}meI m£l m<=I 

\I\=k "-t 1 n <t I 



N 

H^ NtS (x) = H NtS (-x) = Y[ ew(i a -sd Xm ), 5 = +,-. (2.5) 

m=l 

It is readily verified that the new AAOs are also related to the previous ones via 

H- k , s (x) = H N „ k>s (x)H„ N ,s(x), k = l,...,N-l, (2.6) 

and when we set 

#o,5 = 1, (2.7) 

then this relation holds for k = N, too. 

The elliptic kernel function \1/ is now of the form 

%(x, y) = W(x) 1/2 W(y) 1/2 S(x, y). (2.8) 

Here, the weight function is given by 

W(x) = — 1 (2.9) 
C[x)C[— x) 

with C the generalized Harish- Chandra function 

C(x)= TT G(xj-x k -p + ia) ^ 

i</i<iv G ^ + 

The function G(z) = G(r, a + ,a_; z) is the elliptic gamma function reviewed in Ap- 
pendix |Al and the notation 

a=(a+ + a_)/2 (2.11) 
is used throughout Section 2. Also, the special function S is defined by 

3 L k ti g ( x j + p/ 2 ) 

Note that it satisfies 

S(x,y) = S(a(x),r(y)), Vex, r e S N . (2.13) 
Moreover, from the reflection equation flA.61) it follows that 

S(x,y) = S(y,x) = S(-x,-y). (2.14) 
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We are now prepared to recall the kernel identities. They are given by 

(H ltS (x) - H_ l)S (y))*(x, y) = 0, ±l = l,...,N, 6 = +,-. (2.15) 
Equivalently, the 4N commuting AAOs 

Ai iS (x) = W(x)- 1 / 2 H ltS (x)W(x) 1 ' 2 , ±l = l,...,N, 5 = +,-, (2.16) 

satisfy 

(A l>s {x) - A„ l>s {y))S(x, y) = 0, ±l = l,...,N, 5 = +,-. (2.17) 

Using the analytic difference equations (1A.11I) obeyed by the elliptic gamma function and 
the formula flA.20j) relating Rg and sg, it follows that these AAOs have meromorphic 
coefficients. Specifically, one readily obtains the explicit formulas 

A±k,s{x) = 22 II f s ^i x m ~ x n ) 2 Y\_ exj p(Tia-sd Xm ), k = 1, . . . , N, 5 = +,-. 

IC{l,—,N}meI mel 

\I\=k n $I 

(2.18) 

For our purposes, it is crucial that the kernel identities (I2.17P are equivalent to the 
following identities for the functions sg(z): 

IT Sg ~ Xn ~~ ^ J^J s S\ x m ~ Vn + P) 

IC{l,...,N}meI S s( x m-Xn) mg/ Sg(x m - y n ) 

\I\=k n0 ne{l,...,N} 

tt s s (y m -y n + p) tt s s (y m -x n -p) 
ich,N}lh s ^Vm-Vn) H sg(y m -x n ) 

\I\=k n ne{l,...,N} 

This equivalence can be verified by using once more the equations (lA.llj) and (1A.20|) . 
Further details, as well as a proof of (12.191) . can be found in Section 2 of |Rui06] . 

As a preparation for our account of the relativistic Toda regimes we introduce addi- 
tional avatars of the 4N commuting AAOs. They can be defined by 

Af >s (x) = C{Tx)- l A hS {x)C{Tx), I e {±1, . . . , ±N}, 5 e {+, -}, (2.20) 

with C the Harish- Chandra function (I2.10p . so they have meromorphic coefficients as well. 
Alternatively, introducing the elliptic scattering function 

U{x) = C(x)/C(-x), (2.21) 

they are given by 

Af :5 {x) = U{x) ±1/2 H L g{x)U{x)T 1 /\ l = ±l,...,±N, 5 = +,-. (2.22) 
Since we have 

\U(x)\ = l, xeR N , (2.23) 

these operators inherit the formal self-adjointness of the AAOs H^g(x). Note that for 
I = ±N the four sets of commuting operators 

{#,,,}, {A, s }, {A+ s }, {Ai >s }, l = ±l,...,±N, 6 = +,-, (2.24) 

yield the same AAO exp^ias £\ dj). 
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2.2 The hyperbolic case 



In the hyperbolic limit r J, the AAOs H±k,s remain of the same form, but now with 
s$(z) = sinh(7r z / as) , cf. ( 1A.23I) . In addition, the limit (IA.25I) from the elliptic to the 
hyperbolic gamma function implies that the kernel identities (I2.15P hold true if we take 
G{z) to be the hyperbolic gamma function. All other quantities and relations in the 
previous subsection have immediate hyperbolic counterparts as well, so we shall not spell 
them out. 

In the hyperbolic case, however, we are also able to obtain kernel identities relating 
the AAOs 

A ±k>s (x), k = l,...,N, 5 = +,-, (2.25) 

in N variables x — (xi, . . . , xjy) to the following AAOs in N— £ variables y = (yi, . . . , yN-e), 
£ = 0,...,N: 

A T{k ^ s (y), j=0,...,£, 0<k-j<N-£, A 0)S = 1. (2.26) 

In the trigonometric case (which we do not consider), the two hyperbolic periods 
ia + , ici- are replaced by one imaginary period ia and a real period ir/r, and accordingly it 
suffices to consider AAOs Ai, . . . , A N , with the coefficient building blocks 5,5 replaced by 
the sine function. For this case Komori et al. |KNS09] first arrived at the analogs of the 
extra kernel relations for Ai, by using corresponding functional identities. Our reasoning 
below yields relations for arbitrary k, whose trigonometric analogs (with the hyperbolic 
gamma function replaced by the trigonometric one) can be obtained by adapting our 
hyperbolic arguments. 

The relations involve coefficients cf j with £ G N, j G Z, 5 = +, — , given by 

cg )0 = l, 4 = 0, j>£, j<0, (2.27) 

and 

4 = SiMO? - l)p), e s {(i - 3)p), . . . , e s (-(£ - l)p)), j = 0, . . . , £, (2.28) 

where Sj(a\, . . . ,ai) denotes the jth elementary symmetric function of ai, . . . , a^; also, 
here and below we use the abbreviation 



es(z) = exp(nz/a$), 5 = +,-. (2.29) 
Notice that the coefficients are even in p and satisfy 

c e,o = c e,i = 1) c i,j = c \/-ji J = 0) • • • ; ^- (2.30) 
Moreover, it is not hard to verify that the coefficients obey a recurrence relation 

4 + i, = es{jp)4j + e *to - 1 - t)p)4j-i> ( 2 - 31 ) 

and that they are uniquely determined by this recurrence together with the side condi- 
tions ( jzzzp . 

With G(z) denoting the hyperbolic gamma function and A±k t s the hyperbolic version 
of the elliptic AAOs (12.18p . we are now prepared to state and prove the pertinent relations. 
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Theorem 2.1. For £ = 0, 1, . . . , N, let 

^-^nn^^. ^ = i- (2.32) 

AA AJ - G(x m - y n + p/2) 

m=l n=l it/ 

For any k e {1, . . . , AT} and r, 5 G {+, — } we nave 

min(fc,£) 

A Tk;S {xi, x N )S e (x, y) = ^ 4,i^-T(fc-i),5(yi, • • • , yN-e)Se{x, y), (2.33) 

j=o 

waere 

A^fo/i, • • • , J/jv-^) =0, m> N — £, A 0j5 = 1, (2.34) 
and where the coefficients cfj are given by ( 12.271) and ( 12.281) . 
Proof. Since we have 

S £ (x, n) = S*(-a;, -n), (2.35) 

it suffices to show (12.331) for r = +. Our proof proceeds by induction on £. The case £ = 
amounts to the hyperbolic version of (12.171) . so we now assume (I2.33j) for £ > and show 
its validity for + 1. 

To this end we begin by deducing from the asymptotics of the hyperbolic gamma 
function (cf. (1A.32|) and (1A.36|) ) that we have 

lim (f>(x, y N - £ + A)S e (x, yi, ■ ■ ■ , VN-t-u Vn-i + A) = Se +1 (x, y), (2.36) 

A— >oo 



where 

4>(x, z) = exp I y~)(x m - z)) ■ (2-37) 

\ m=l 



In order to exploit this limit, we note 

(j}(x,y N - e )A k>s {x) = e s (-kp)A k ,s(x)<f)(x,yN-t)- (2.38) 

Furthermore, we split the A AOs A_ m ^s(y) into two parts, depending on whether the index 
set I contains N — £ or not: 

A- m , s (y)= £ (•••)+ £ (•••)• (2 ' 39) 

IC{l,...,N-£-l} I=JU{N-£} 

\I\=m Jc{l,...,N-l-l} 
\J\=m-l 

Denoting the first and second sum by B m> s(y) and C mt s{y), respectively, we then observe 
that 

4>(x,y N -e)B m: s(y) = B mi s(y)(/>(x,y N -e), 

<j>(x,y N -t)C m ,s{y) = e 5 (-Np)C m:5 (y)<i)(x,y N - e ). (2.40) 

Next, we multiply both sides of ( 12.331) by the function es(kp)(fi( y x,yN-e) and use the 
commutation relations (I2.38j) and (12.401) . We then take yN-i Un-e + A, and use the 
readily verified limits 

lim B mS {y u . . . , y N -e-i, Vn-i + A) = e s (-mp)A_ mS (y 1 , . . . , y N -£-i), (2.41) 

A— >oo 
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lim C m>s {yi, Vn-l + A) = e s ((N - £ - m)p) exp(ia- S d y J 

A— >oo 

xA^y(i/i, ...,|/jv-^-i), (2.42) 

together with the limit (12.361) . If we now take j — > j — 1 in the sum coming from Ck-j t s(y), 
then it becomes clear that the coefficients of the AAOs A^fk-j),s on the right-hand side 
satisfy (ESQ), so that they are given by (f2T28l) with t->l + l. □ 

We proceed to detail three specializations of Theorem 12.11 For the first we fix k = 1, 
but impose no restrictions on £. From (I2.28P we obtain 

4 = 1, 4 tl = s s (£p)/s s (p), (2.43) 

yielding the following special cases. 

Corollary 2.2. For £ = 0, . . . , N , we have 

ss(£p) 

(A± ltS (xi, ...,x N )- A TltS (yi, . . . , y N -i))S e (x, y) = — f-r-S t (x, y). (2.44) 

We note that this corollary is the hyperbolic analog of Statement (1) in Theorem 2.2 
of |KNS09j . 

Next, we require £ — 1, but do not restrict k. From (12.28P we have c\ = cf 1 = 1, and 
hence the following specialization results. 

Corollary 2.3. For k = 1, . . . , N , we have 



= (AfkAvii • • • ' 2/Jv-i) + • • • > Vn-i))Sx(x, y). (2.45) 

Finally, we choose £ = N, recalling Sn = 1. 
Corollary 2.4. The following functional identities hold true: 

E-tt s s (x m - x n ± p) 5 . . 

I(Z{X,...,N} mel 
\I\=k n P 

These identities were obtained before in the proof of Lemma A. 5 in |Ru i95j. 



2.3 The periodic Toda case 

As explained in [Rui94j . the periodic Toda AAOs ( ll.20p can be obtained as limits of the 
elliptic AAOs ( II. 9p . In this section we shall in particular recover this result as a corollary 
of somewhat more general limit formulas, detailed in Lemma 2.5. More precisely, these 
formulas are primarily derived to obtain Toda kernel functions, but they can also be 
used to show that the AN commuting AAOs H± k ^(x) given by (12.ip - (j2.3p give rise to 
AN commuting periodic Toda counterparts, denoted by the same symbols. As it shall 
transpire, however, the AAOs A±k,s ( I2.18P have no sensible limits (the N = 2 case being 
a curious exception). 
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This makes it all the more surprising that the elliptic functional identities (I2.19P 



corresponding to the relation between S(x,y) and A^s (as expressed in ( 12. 17ft ) do have 
Toda limits. Once this limit is obtained, we can easily identify kernel functions S' ± (x,?/) 
for the periodic Toda system. But here these functions correspond to AAOs that are the 
periodic Toda counterparts of the elliptic operators Af 5 given by (I2.22p . Indeed, as in the 
elliptic case, they are the similarity transforms of the Toda AAOs His with a function 
U(x). For a suitable choice of parameters this Toda [/-function is unitary. Moreover, this 
parameter choice entails that all of the Toda AAOs H^s and Af s are formally self-adjoint. 
Turning to the details, our starting point consists in making the substitutions 



and 



X n -»■ X n - — , Vn-^Vn - "TT i u — 1, . . . , N, (2.47) 

JMr Nr 

p^p+ih> (2 - 48) 

in the elliptic quantities occurring in Subsection 2.1. Next we consider the limit r 
0. (Recall the real elliptic period 2u> is parametrized as n/r.) It seems intractable to 
control this limit for the quantities expressed in elliptic gamma functions, and in fact it 
appears likely that none of them can be renormalized so that this limit exists. Rather, we 
concentrate on the functional equations (I2.19p . which only involve the functions s±(z). 

We first note that the product representation ( 1A.21I) for the function sg(z) contains an 
exponential factor e$(— rz 2 /ir). Rather than directly performing the substitutions ( 12. 47ft 
and (12.481) in (12.19]) . we may and shall eliminate these factors from the start. The point 
is that the identity ( 12 . 1 9[) still holds true if we switch from s$ to the function 



s s (z) := e s (rz 2 /n)s s (z) 

oo 

smh{nz/a$) J~| 



a 5 . , , , , tt (1 - e s (2z - 2vrZ/r)) (z -»■ -z) (2.49) 



7T 



1=1 (l-e s (-2nl/r)Y 



Indeed, a straightforward computation shows that the exponential factors combine to 
yield the same overall factor in the left-hand and right-hand side. 

Accordingly, we substitute (I2.47[) and (I2.48[) in (12. 19[) with ss replaced by sg, and 
proceed to study the asymptotic behavior as r — > of the resulting identity. To this end 
we focus on the factors s$ (y m — x n — p)/ss(y m — x n ). Their asymptotics is given by the 
following lemma, which involves an auxiliary Toda building block 

t s (z) = l-e s (2z + 2p). (2.50) 

Lemma 2.5. Let m, n = 1, . . . , N, and m ^ n. Then we have, as r — > 0, 

ss (y m ~ x n ~ P + -jfc{n - m - 1)) 
ss (y m - %n + jfc(n - m)) 

( m-n ( 7T \\ J t s (x m+ i -y m ), n = m + 1 (mod N), , , 

e& \\m- n\ V Nr) I I 1, otherwise. ^' bL) 



Moreover, 

~"(y™- x m-p-ib) „ e*(P_+jk) r _, . (2 . 52) 
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Proof. We shall infer the statement from the product representation (I2.49P for the func- 
tion s$. Dealing first with the sinh-prefactor, we readily find 

sinh f g (y m - x n - p + jfc(n - m - 1)) 

sinn t s (y™ ~ Xn + jk( n ~ m )) 

„ JpZ± ( p+ lS] fts(x m+1 -y m ), n = m + l, 
\\m — n\ v Nr/ J y 1, otherwise, 

as r — > 0. Next, we consider the terms arising from the infinite product in (12. 49ft . Since 
|n — m\ < N and I > 1, we have 

1 - e s (2(y m - x n - p) + (*=%=± - I)) 
lim V , — 7 9 r / * — = 1, 2.54 

and 

1 - e, (-2(y m -x n -p) + * (-^ ~ 0) 

r^O 1 - es {-2{y m - Xn ) + f (-^ - 0) 

ts(x 1 -y N ), m = N,n = 1,1 = 1, , . 
1, otherwise. 

Putting the pieces together, we arrive at the statement for m ^ n. The remaining case 
m = n now follows easily, noting all factors in the infinite product in (I2.49P then converge 
to one. □ 

It is clear from (I2.49[) that Sg(z) is an odd function. As a corollary of (j2.51j) . we thus 
have 



ss 



{x m -y n + p+ j^{n - m + 1)) 



ss (x m -y n + j^{n - m)) 



n-m / tt \\ j t s (x m -y m ^), n = m - 1 (mod N), 
65 I 1 P + iV>/ J 1 1, otherwise. l ^ bbj 



\n — m 



In addition, if we replace y m by x m in (12.511) . and x m by y m in (12.561) . then we clearly 
obtain the contribution to the asymptotics due to the remaining factors in (I2.19p . 

With these asymptotic formulas at our disposal, we are in the position to obtain the 
following Toda counterpart of the elliptic functional identities (I2.19p . 

Lemma 2.6. For k = 1, . . . , N , we have 

ts(.Xm Vm— l) 



/■■ (i v» ,„-/ mei ts ( ym Xm P ^ 

\I\=k m+l<£I 



/_!! \] m &i t& ( ym Xm p } 

\I\=k m-l$I 



withtg(z) given by (I2.50p . 
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Proof. Just as for (|2.19p . the special case k = N of (I2.57P is obvious. Fixing k < N, 
we substitute (12 .47ft and (I2.48j) in the identity ( 12. 19ft with ss — > sj, and then exploit the 
above asymptotic formulas in the following way. We focus on a term in the sum on the 
left associated with a fixed subset I. First, consider a pair of indices m 6 / and n I. 
This gives rise to a product of ratios 

5g(x m -x n -p) s s (x m -y n + p) 

■5<5^-£m ^nj SS\X m Vn) 

With the substitutions in place, we can now use (12.511) for the first and (12.561) for the 
second ratio to deduce 

{ts(x m +i- x m ), n = m + 1 (mod N), 
ts(x m -y m -i), n = m - 1 (mod N), (2.59) 
1, otherwise. 

The remaining pairs of indices m, n £ / yield the product 



n Ss{ -t ! n V } = n t - n t — ( 2 - 6 °) 



m>n 



where we have introduced the ratios 



m s s{x m Vm~^~P) fn 

m = ~~H — Z — S~> ' ' 

S8\X m Vra) 



and the product of ratios 



rj, s s {x m -y n + p)s s (x n -y m + p) 

T m n = — ^ = 7 — ( = \ ■ (2-62) 

SS\X m Vn) SsyX n y-m) 

In view of f !2.52j) we have 

lime 5 \ -p - — J T m = l/t s {y m -x m - p). (2.63) 

r->0 V ISr/ 

Moreover, using oddness of §s in the second ratio of T mn , we can use (I2.51I) for the second 
and (I2.56I) for the first ratio to get 

{ts{x m -y m -i), n = m-l, 

ts(x!-y N ), m = N,n = l, (2.64) 
1, otherwise. 

A moment's thought now shows that when we multiply the left-hand side of (12.191) 
with sg — > ss by the renormalizing factor 

then its r — > limit yields the left-hand side of (I2.57p . Proceeding in the same way for 

the right-hand side of ( 12.191) with ss —> §s, we then deduce ( 12.571) . □ 
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At this point we invoke the modified hyperbolic gamma functions Gr and Gl, cf. Ap- 
pendix [A] Indeed, the identities (I2.57p . combined with the difference equations ( 1A.34I) 
and ( 1A.35I) satisfied by Gr and Gl, resp., can now be used to obtain kernel functions for 
the periodic Toda system. Specifically, it is readily deduced from these formulas that we 
have 

II t f Xm - ym 2 1 \ = AHx,y)- 1 ]\expteia^a x JA ± {x,v), (2-66) 

and 

n y^ 1 "^ = A±( X , y)- 1 n ex P (± M „,^)^(,, „), ( 2 . 6 7) 

where we have introduced the auxiliary functions 

N 



11 G L (y m - x m -ia) 



m=l 

= ft C t (». -*, + *.) 
~y (jRyVm ~ Xm+l +ia- p) 

These functions are basically the kernel functions we need. To detail this, we should first 
introduce various additional Toda quantities. 

To begin with, we have thus far retained the coupling parameter p of the elliptic 
and hyperbolic regimes, since this yields the simplest auxiliary quantities. At this stage, 
however, we need to switch to parameters that are more appropriate for the Toda regimes. 
First, we introduce four Toda interaction functions 

Tf(z) = \ + e B (2z±ia- S + 2-q), 5 = +,-, 77 e K. (2.70) 

They are obtained from the four functions tg(z) and tg(z — ia_s) when p is replaced by 
ia+rj. (The function T+(z) amounts to the interaction function (ll.2ip of the Introduction, 
cf. (jl.23p .) The real parameter rj plays the role of coupling constant. We also introduce 
the Toda ^/-function 

n 1 

m=l 

It has the unitarity property 

\U(x)\ = l, xeR N , rjeR, (2.72) 

and it is invariant under cyclic permutations of x%, . . . , Xn- 
Next, we define a set of AN formally self-adjoint AAOs 

At,si x ) = II T ^( x m+i-x m )Y[exp{-ia_ s d Xm ), (2.73) 

\I\=k mel mel 
m+lgl 

At kfS (x) = J2 II T ^ x ™ - x ™~^ II e*P(ia-sd Xm ), (2.74) 

\I\=k m£l m£l 
m—lfil 
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and a second set of 4N formally self-adjoint AAOs 



A,s( x "> = Yl n T <5 ( Xm ~ Xm -^ II exp(-m__ 5 9 Im ), (2.75) 



\I\=k mel mel 



A -k,s( x ) = Yl II T & ( x m+i-x m )W^{ia. 5 d Xm ). (2.76) 

\I\=k m&I m&I 

It is easy to verify that these definitions entail the relations 

N 

A T _ N , 5 {x) = A^six)- 1 = J] exp(za_ 5 ^J, (2.77) 

m=l 

^(x) = A T N _ k!S (x)AL N;S (x), k = l,...,N-l, (2.78) 

where t,5 = +,—. Moreover, from the difference equations (IA.350 obeyed by Gl we 
obtain 

Al 5 (x) = U(x)Al 5 (x)U(x)-\ ±1 = 1,..., N, 5 = +,-. (2.79) 
Finally, we introduce a third set of 4iV formally self-adjoint AAOs, namely, 

H liS (x) = Uixf^Al^Uix)- 1 ' 2 = U{x)-^ 2 Al s {x)U{xf' 2 , ±l = l,...,N, 5 = +,-. 

(2.80) 

These operators are the quantum counterparts of the classical Hamiltonians given by ( 11. II) . 
( 11.151) and ( 11.161) . with modular invariance and formal self-adjointness taken into account. 
Explicitly, letting k = 1, . . . , N, and 5 = +, — , they read 

*m*)=z) n r/(vi-^) i/2 n^ A " n ^^-x^y 2 , ^ 

\I\=k m£l m&I m&I 

H- h j{x) = Y, II T-{x m+1 -x m ) l ' 2 J[e ia -^ J] Tf(x m - x m ^f' 2 , (2.82) 

as is readily checked by using (1A.35|) once again. It is also clear that they satisfy (12.771) 
and (I2.78P with A T replaced by H. Furthermore, like the Toda AAOs A[ s , they are 
invariant under cyclic permutations. 

Last but not least, the 4iV operators H^g mutually commute, so that the sets of op- 
erators {*4/ 5 } and {.A^} consist of mutually commuting operators as well. As explained 
in the Introduction, this assertion follows from [Rui90j . where a direct proof of commu- 
tativity can be found. 

On the other hand, taking commutativity of the elliptic Hamiltonians H^g in Subsec- 
tion 2.1 for granted, the commutativity of the Toda Hamiltonians H^g also follows from 
the latter being limits of the former. A quick way to check these limits within the present 
context is as follows. First, push the third product in (12. ip through the shifts. Then re- 
place Sg by Sg, which amounts to a multiplicative renormalization. Now substitute (12.471) 
and (I2.48P and use (I2.5ip and (I2.56P (with y replaced by x) to see that the r — > limit 
yields flZEED - 

The reader who has verified these steps will easily see why this procedure fails for the 
AAOs Ak t s(x) with k < N, unless N = 2. The point is that different subsets / in (12.181) 
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give rise to different powers of the factor es(iT/Nr), so that no nontrivial r — > limit can 
be obtained by a multiplicative renormalization. 

Possibly, the substitutions ( 12.47}) and ( 12. 48ft in the elliptic AJ S , along with a suitable 
renormalization and similarity transformation, yield AAOs that converge to the periodic 
Toda AJ S for r — > 0. At any rate, the analogous substitutions in the hyperbolic Aj s do 
yield the nonperiodic Toda AJ S as limits, cf. the end of the next subsection. We have 
anticipated this state of affairs in the notation we have adopted above. 

We are now prepared for the main result of this subsection. 

Theorem 2.7. Let I G {±1, • • • , ±N}, and 5 G {+, — }. We have kernel function identi- 
ties 

(A+ s (x)-At l!S (y))S + (x,y) = 0, (2.83) 
(Ai/x) - AZ llf {y)) S-(x, y) = 0, (2.84) 
(H l>s (x) - H- l)S {y)) UixT^Uiyr^S+ix, y) = 0, (2.85) 
(Hu(x) - H-M) U{x) l / 2 U{yfl 2 S-{x, y) = 0, (2.86) 



where 



AJ i G L (y m -x m + ia/2 + r l /2) 1 ; 



m= 

A' 



^Ji G R {y m - x m+1 + ia/2 - r]/2) 

and U(x) is given by (12. 71 j) . Furthermore, the identities (12. 83p - (12.861) still hold when the 
functions S ± (x,y) are replaced by S ± (y,x) or by S ± (a(x), y), with a any cyclic permuta- 
tion. 

Proof. Combining the functional equations (12.571) with (12.661) and (12.671) . we obtain 
Y 1 1 t si x m+i ~ x m) Yl exp(-ia_ s d Xm )A + (x, y) 

IC{1,...,N} mel mGl 
\I\=k m+l£I 

= Yl II hiVm-Vm-i) Y[ew(ia-sd ym )A + (x,y). (2.89) 

IC{1,....N} m£l m£l 
\I\=k m-lfl 

If we now replace p by ia + r\ and then shift y\, . . . ,y^ by 3ia/2 + r//2, then we obtain 
(12.831) for I = k. For I = — k, we can use (12. 77p (12. 78j) . together with the relation 

AZ N /x)S + (x, y) = A^ 8 (y)S + (x, y), (2.90) 

to complete the proof of (12.831) . 
Likewise, from the identities 

^ ] [ t 5 (x m+1 - x m ) Y\ exp{ia-6d Xm )A~~(x,y) 

IC{1....,N} m&I mGl 
\I\=k m+10 

= II hiVm-Vm-i) Ylexp(-ia-5dy m )A~(x,y), (2.91) 

IC{1,...,N} m£l mel 
\I\=k m-lil 
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we obtain f |2.84[) with I = —k upon replacing p by — ia + rj and shifting y%, . . . , y^ by 
-3za/2 + rj/2. Then 11233]) with / = jfe follows as before. 

Recalling (ESQ]), we now obtain (l2T8fl -( l2~86]) from (l2~83]) -( l2T84l) . Also, the last state- 
ment follows from invariance of the AAOs and U (x) under cyclic permutations. □ 

In the Introduction we have already pointed out that from a given kernel function we 
can obtain an infinity of other ones, cf. the paragraph containing f 1 1 . 2 9 j) . The kernel func- 
tions in Theorem 12.71 however, are not related to each other by a coordinate translation 
or multiplication by a factor (jl.29p . 

It is worth pointing out that limits of translations lead to elementary kernel functions, 
due to the simple asymptotics of Gr(z) and Gl(z) for $l(z) — > ±oo, cf. ( ]A.36j) and( 1A.37l) . 
To be specific, consider the substitution 

y m ^y m -ia/2-r)/2 + A, m = l,...,N, (2.92) 

in S + (x,y). For A — > oo, the dominant asymptotics is given by 

N ( in \ 
exp(2iiVx) TT ex P ~^-(fm - x m + A) 2 . (2.93) 



Thus, if we multiply by a factor 



(2i7rA N 
V(x m -y m ) ) , (2.94) 



m=l 



(which is of the form (I1.29P ). then we can take A — » oo and conclude that the function 



N 

1% 



k r( x ,v) = ex P y^( x m ~ V m f , (2.95) 



m=l 



is a kernel function for the AAOs A.t s . Likewise, the left asymptotics yields an elementary 



kernel function 

N 

ITX 



Kt(x,y) = exp V(a; m - y m ^) 2 . (2.96) 

» a+a_ / — ' ' 

m=l 



Notice that the latter results from by a cyclic permutation. 
In the same way we obtain elementary kernel functions 



N \ / . N 

— 111 \ , , o\ , / — Z7T \ , 2 



K R ( x ,y) = ex p ( ^-X^^ m_2/m ) 2 ) > K L( x >y) = ex P ( 7-7- $^( Xm ~ Vm - x ' 

(2.97) 



a_i_a_ ^ — ' I \ a + a_ 

m=l / \ 1 m=l 



for the AAOs A l s . 



2.4 The nonperiodic Toda case 

In this subsection we first deduce kernel functions for the nonperiodic Toda system by 
limit transitions from the periodic case. We then take further limits to obtain kernel 
identities that relate AAOs whose number of variables differs by one. Finally, we discuss 
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the direct limit from the hyperbolic quantities to their nonperiodic Toda counterparts, 
which is quite easily understood for the AAOs AJ s (x) as well. As a bonus, we obtain the 
nonperiodic Toda [/-function as a limit of the hyperbolic one. Throughout this subsection 
we use the same symbols for the nonperiodic quantities as for their periodic counterparts. 
We start by observing that if we perform the substitutions 

x n -tx n + nA, n = l,...,N, rj -»■ rj — A, (2.98) 

in the AAOs Aj s (x), then only the factor T$(xi — xn) is affected, and it converges to 1 
in the limit A — > oo, cf. (12.701) . Since we retain the notation Aj s (x) for the AAOs thus 
obtained, the equations ( I2.73l) -( j2.76l) still hold true, but now with the nonperiodic Toda 
convention xq = —xn+i = oo. Likewise, on account of the G^-asymptotics ( 1A.36I) . we 
have the relation ( 12. 79ft between Af s (x) and Af s (x), but now with 

N-l 

U{x) = TT — -. (2.99) 

m=l 

Finally, we have nonperiodic Hamiltonians Hi^(x) given by (12.80p - (12.82j) . It should be 
noted that the nonperiodic AAOs are no longer invariant under cyclic permutations of 
the coordinates xi, . . . ,xn- 

Next we discuss the effect of substituting first 

Vn -> Vn -v/2 + nA, n = 1, . . . , N, (2.100) 

and then (12.981) . on the kernel functions S ± (x,y) given by f!2.8Tj) (I2.88P . Clearly, the 
resulting functions depend on A only via the factors Gr{ijn — x\ ia/2 — r] + NA). It 
follows from ( 1A.361) that if we take A — > oo, then these factors converge to 1. Shifting 
next y n back by r]/2, we wind up with kernel functions 

c+ rr v) = \ TT GR{ym-x m ^-ia/2--n/2) 

1 ' y) G L (y N -x N + ia/2 + rj/2) 11 G L (y m -x m + ia/2 + V /2) ' 1 ' ] 

Q-f \ n ( ■ m a /o\ TT G L (y m ~x m -ia/2 + r]/2) 

S (x,y) = G L (y N - x N - la 2 + i] 2) I I — — — — — . 2.102 

^ G R {y m - x m+1 + ia/2 - rj/2) 

We thus arrive at the following analog of Theorem 12.71 for the nonperiodic Toda system. 

Theorem 2.8. With the convention 

xo = yo = oo, x N+ i = y N+1 = -oo, (2.103) 

in force, the identities (12T83|) - d2~86l) hold true for S ± given by (12~T0T|) - ([2302]) and U 
by (12T991 . 

Just as in the periodic Toda case, the substitution (I2.92|) and subsequent limit A — > oo 
lead to the elementary kernel identities 

(Al s (x)-AL l>s (y))K- R (x,y)=0, ±l = l,...,N, S,r = +,-, (2.104) 
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where 

A^(x,y)=exp( — f> m -y m ) 2 ) . (2.105) 

\ a + a ~ m =i / 

By contrast, in this case the limit A — > — oo does not yield the kernel functions 
defined by (12. 96ft and (12.97!) . due to the 'missing' G^-factor. In fact, since the nonperiodic 
AAOs are not invariant under cyclic permutations, one should not expect that (12.1041) 
also holds for K T L . 

We proceed to a more significant difference between the periodic and nonperiodic 
Toda systems: from Theorem 12.81 we are able to deduce kernel identities that relate 
the AAOs Aj s in A" variables x = (xi, . . . , xjy) to AAOs Aj, s in N — 1 variables y = 
(y 1 , . . .,y N -i). 

Corollary 2.9. Letting k G {1, . . . , N}, and 5, a, r e {+, — } ; we have 

Al k j{xi, . . . ,x N )Sl{x,y) 

= {^-ak,s(y^ • ' • ' ftf-i) + A -*(k-i)Ayi> • • • ' Vn-i)) SI(x, y), (2.106) 



where 
and 



A T ±NiS (yi, . . .,y N -i) = 0, Als = 1, (2.107) 



iV-l 



St(x,y) = TT G * {V r ~ Xm+1 " "I 2 " # , (2-108) 
' AJ i G L ( 2/m -x m + za/2 + r / /2) 1 J 



m= 
JV-1 



S-(x v) - TT ^fo" - g " ~ ifl /2 + q/2) f , 

Proo/. Clearly, fTA~36|) entails that if we substitute y N ^ y N - A in (12.1011) and (12.1021) 
and let A — > oo, then we obtain S^(x,y). In order to determine the same limit for the 
AAO 

AL ak>s ( yi ,...,y N ), (2.110) 

we split the sum in the pertinent formula among (I2.73I) - (I2.76I) into a sum over subsets I 
containing the index A^ and a second sum over I not containing N. Now there are two 
cases to consider. 

(1) ( a — t) If k — N, then the second sum is empty, whereas for k < N it equals 
A T _ Tk s (y\, . . . , Vn-i), cf. (I2.74p - (l2.75p . This AAO is independent of yN, so it is invariant 
under the limit at issue. On the other hand, if we perform the substitution in the first 
sum and let A — > oo, then we obtain the AAO 

A -r(k-i),d(y^ • • • ' ftf-O exp(Tia_ s dy N ). (2. HI) 

Since Sl(x,y) is independent of y^, the AAO exp(rm_5(9 yjv ) acts as the identity on 
Sl(x,y), and so we arrive at ( 12.1061) . 

(2) (a = — r) Arguing as before, we see from (12.731) and (I2.76P that in this case the 
roles of the first and second sum are reversed, i.e., they now yield the first and second 
AAO on the right-hand side of (12.1061) . respectively. □ 
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It is worth pointing out that we cannot repeat the above limit procedure for S{(x,y) 
and the variable yjv-i, so as to obtain a kernel function Sl{x 1 y) that relates AAOs in N 
and N — 2 variables. Indeed, the asymptotics ( 1A.37I) of Gr and Gl implies that if we were 
to renormalize S"[(x,y) so as to obtain a finite limit, then the resulting function would 
contain an exponential factor that still depends on Vn-i- Likewise, it seems impossible to 
obtain analogs of S^(x,y) for the periodic Toda case via a limit of the kernel functions 
S ± (x,y) given by <^M>-<^M). 

To conclude this subsection, we discuss an alternative way to arrive at the above 
nonperiodic Toda results, namely, via the hyperbolic quantities. The pertinent limits are 
far more easily controled than the elliptic to periodic Toda limits, and yield additional 
insights. First, we substitute 

x n -> x n - nA, y n ^y n -nA, p -> p + A, (2.112) 

in the hyperbolic version of the functional equations (12.191) (obtained upon replacing ss(z) 
by smh(jrz/as))- The asymptotic behavior is then given by 



sinh f s (Vm ~ x m - p - A) es ( p + A) 



sinh^(2/ m -x m ) U(y m - x m - p)' 



A^oo, (2.113) 



and by ( 12.531) with n/Nr replaced by A (recall tg is defined by (12.501) ). 

Next, we use these formulas in the same way as in Subsection 2.3 to arrive at the 
functional equations (12. 57ft with the nonperiodic Toda convention (12. 1Q3[) in effect. More 
specifically, in the proof of Lemma [2761 we need only replace Ss(z) by sinh(7Tz/as), and the 
r — > limit by the A — > oo limit. 

The developments below Lemma 12.61 can now be followed verbatim, with the con- 
vention ( 12.103!) ensuring that we get the same quantities as obtained via limits of the 
periodic Toda regime. But now we can also control the limits of the hyperbolic AAOs 
and U-f unction. 

Indeed, using the difference equations (1A.24I) obeyed by the hyperbolic gamma func- 
tion, we can readily calculate the hyperbolic AAOs Aj s explicitly from (I2.10p . (I2.18P and 
flZZDD - This yields 

■^Xs( x ) = ^2 IT U-( x m - x n ) 2 f s ,+ (x m -x n - ias) 2 Y\_ ex P(-i a -sd Xm ), (2.114) 

\I\=k mel,n^l mel 
n>m 

<At k j(x) = ] [ fs,+ (x m - x n ) 2 f S -(x m - x n + ia^ s ) 2 Y\_ ex P( ia -s d xrn)i (2.115) 

\I\=k mel ,n^I mel 

Afc.iO'O = ^2 FI h-( x m ~ x n ) 2 f s , + {x m - x n - ia^ s ) 2 JJ exp(-ia_ s d Xm ), (2.116) 

\I\=kmeI,n£I mel 

•AZ k ^{x) = ^2 FI fs,+ ( x m- x n ) 2 f 5 -(x m ~ x n + ia^ 5 ) 2 Y\_exp(ia_ 5 d Xrn ), (2.117) 

\I\=k m£l,n(£I 
n>m 

where k = 1, . . . , N, 5 = +, —, and 



\I\=kmeI,n$I mel 
n>m 



smh(Tt(z ± p)/as) 



fsA*) = -A i \ • (2-118) 

smh(7rz/ ag) 
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Now we substitute f 1 2 . 1 1 2 j) and take A to oo. Denoting the limits by the same symbols, 
this yields 

Atsi x ) = ^ 1 [ t s (x m+1 - x m ) Y\ exp(-ia- S d Xm ), (2.119) 



\I\=k mel mel 



^-kA X )=^2 H t5 ( Xm ~ Xm ~^ II eyj P( ia -S d x m ) , (2.120) 



\I\=k mel mel 



\,s( x ) = n U{x m - x m ~i - ia-s) Y\_ ex P(- ia -sd Xm ), (2.121) 
\l\=k mel mel 

m—l^I 

AZ k ,6( x ) = ^ ] I t s (x m+1 -x m - ia_ s ) Y[ exp(ia- S d Xm ), (2.122) 

\I\=k mel mel 

with t$(z) given by (I2.50p and the convention (I2.1Q3[) in effect. 
If we now set 

p = ia + rj, (2.123) 

then the AAOs (I2.1190 - (l2.122j) turn into the nonperiodic Toda AAOs already obtained 
via the periodic regime. In the same way, the hyperbolic Hamiltonians H±k.s converge to 
their nonperiodic Toda counterparts, whereas the hyperbolic AAOs A±k,s have no sensible 
limit, a feature shared by the hyperbolic weight function W(x). 

The hyperbolic scattering function U(x), however, does have a limit, provided a suit- 
able renormalization is performed. To be specific, when we substitute (12. 1 1 2j) in U(x) 
(given by (I2.2ip and (I2.10p with G the hyperbolic gamma function), then we obtain via 
the G-asymptotics (which can be gleaned from ( 1A.31|) . (1A.32|) and (1A.36j) ) 

JV-l 

lim N N{N - l ^ 2 U(x) = TT —. -, (2.124) 

a^oo J-j. G L (x m+ i - x m + p - ia) 

m=\ v ' ' 

where M is the renormalizing constant 

(171 \ 
{(p-ia + Kf + a 2 ) . (2.125) 
a + a_ / 

If we now again replace p by ia + rj, then the limit function on the right-hand side of 
(I2.124p turns into the nonperiodic Toda [/-function given by (I2.99p . as announced. 



2.5 The dual nonperiodic Toda case 

As recalled in the Introduction, the hyperbolic relativistic Calogero-Moser system is self- 
dual. In the limit leading to the nonperiodic Toda dynamics, which we have just discussed, 
this self-duality property is not preserved. But the dual Toda quantities can be obtained 
by a similar, but simpler scaling limit. Turning to the details, we substitute 

p^p + A, (2.126) 

and study the limit A — > oo. As explained in the paragraph containing ( 11.28}) . this 
should be done for the hyperbolic quantities expressed in the 'spectral variables' v±, . . . ,vn 
instead of the 'geometric variables' x±, . . . , xn- 
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For the resulting dual hyperbolic AAOs 

a f\ ST TT sinhOKn - v n T p)/a s ) tt a , , 010 ^ 

|J|=fcmeJ toG/ 

we get finite limits by a renormalization. Specifically, we readily deduce 
lim e 5 (-k(N - k)A)A± k>s (v) = 

A-j-oo 

e 5 (k(N - k)p) II o 5{ 7} V ( m \ Vn \ II eMTiasd Vm ). (2.128) 

\I\=km£l v me/ 

Here and from now on we use the notation 

s s (z) = sinh(7T2;/a 5 ), 5 = +, (2.129) 

which already occurred in Subsection 2.2. 
Next, we introduce the product 

p(v)= n ex p(^v^-^) 2 V (2- i3 °) 

l<j<k<N ^ + ' 

Using 

P{vY x eMTia- S d Vm )P(v) = J] e s (±(v m - v n ) - ia_ 5 /2), (2.131) 

it becomes clear that, up to a /c-dependent multiplicative constant, the limit AAOs on the 
right-hand side of (I2.128P are the similarity transforms under P(v) of dual AAOs defined 
by 

A ±k ,s(v) = ( T z)^- fe ) II o 7 S II exp(Tia_ a a Um ). (2.132) 

\I\=km£l rnel 

Here we have k = 1, . . . , N , 5 = +, — , and the phase choice will be clarified shortly. 

Consider now the kernel function S(v,w) given by (12.121) . with G the hyperbolic 
gamma function. Shifting Vj by p/2, substituting (12.1261) . and letting A — > oo, it follows 
from the G-asymptotics that the dominant term is a product of 

JV 

S{v,w)= 11 G( Vj -w k ), (2.133) 
j,k=i 

and A-dependent quadratic exponentials. Likewise, when we shift Vj by —p/2, substitute 
(I2.126p . and let A — > oo, we obtain 1/S(v,w) times A-dependent exponentials. Remark- 
ably, when we omit the exponentials, we wind up with kernel functions for the dual AAOs, 
as shown in the next theorem. 

Theorem 2.10. Letting I G {±1, . . . , ±iV}, 5 G {+, — } and a G {±1}, we have the dual 
kernel function identities 

(Au(v) - A_u(w))S(v, wY = 0. (2.134) 
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Proof. From f |2.132[) we see that 

A l>s (-v) = A_ lfi (v), (2.135) 

whereas the reflection equation for the G-function entails 

S(-v, -w) = 1/S{v, w). (2.136) 

Therefore, it suffices to show (12. 134f) for a — 1. 
Now it also follows from (12.1321) that 

A- kt s{v)=A N - k)S (v)A- Ni s{v), k = l,...,N. (2.137) 

Hence we need only show (12.1341) for a = 1 and I — k > 0. To this end, we invoke the 
difference equations (IA.24I) for the hyperbolic gamma function to obtain 

1 N 1 

5(1;, wr eM-ia-sd v JS(v, w) = ^ ]J ( _. fl) , (2.138) 

71=1 y 

1 N 1 

S(v,wr"eMia-sd w JS(v, W ) = ^ \{ _ - . (2.139) 

n=l ' 

From this we deduce that we are done if we can prove the functional identities 

v n - n - 

\I\=k nil ne{l,...,JV} 

= e — v n — t — ~ — (2.140) 

r^ri *n ti S s{ w m ~ W n ) ~~ S 6 {v n - W m ) 

\I\=k nil ne{l,...,N} 

In order to show that these identities are valid, we substitute (12.1261) and x = v,y = w, 
in (the hyperbolic version of) (12.191) and let A — > oo. Then we obtain from equality of 
the dominant asymptotics the identities 

y> TT e S (v n - V m ) -i-r e S {v m ~ W n ) 

rr-n An A S s( V n ~ v m) S S (v m - W n ) 

IC{1....,N} mel m&I 

\I\=h n ne{l,...,N} 

TT e 5 (w m - W n ) yr e S (v n - W m ) 

r ss(w m — w n ) s$(v n — w m ) 

IC{1 ,N}meI v ' m€l ° v n m ' 

\I\=k n<£I ne{l,...,JV} 

Now it is not hard to see that (12.1411) is equivalent to (I2.140p . Indeed, consider the product 
of all e^-factors for a given /. Taking / equal to {1, ... , k}, we get on both sides a factor 



e«5 



' N \ 

*5>i-«>i)) ■ (2-142) 



Since this factor is permutation invariant, it does not depend on /. Hence we can cancel 
all exponentials and obtain (12. 140j) . □ 
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When we make the substitution (I2.126[) in the dual hyperbolic AAOs Aj s (v) (given 
by (I2.114p -f l2.118p ). then the power of e^A) in the dominant asymptotics of each term 
in the sum depends on I. Hence we cannot renormalize these operators so as to obtain 
nontrivial finite limits for A — > oo. By contrast, for the dual hyperbolic Hamiltonians 
H±k,s(v) we obtain 

lim e s (-k(N - k)A)H ±ki5 (v) = e s (k(N - k)p) 

es(T(v m - v n )) \ 1/2 yj . , yj f e 5 (±{v m - v n )) \ 1/2 

me/ ma " ' 

(2.143) 



A— >oo 



En 

\I\=km£l 
nil 



2s 5 {T{v m - V n )) 



Pushing the exponentials through the shifts, we infer that these limits are constant mul- 
tiples of dual Toda Hamiltonians formally given by 



h**m = e n ( 

\ T\ — U rnC T \ 



1 



2s S (T(Vm ~ V n )) 



1/2 



Y[ exp(Tia- S d Vm ) J J 



me/ 



me/ 
n0 



1 



2s S (±(v m - V n )) 



1/2 



\i\=kmei 

(2.144) 

Due to the square root ambiguity, the phases of the individual terms in the sum are 
not well defined. To remedy this, we first note that the relevant Hilbert space is L 2 (G, dv), 
where the dual Toda configuration space is defined by 



G = {v eR N \v N < ■■■ < Vl }. 
Now we fix the phase ambiguities by defining 



H 



±k,S{V) 



En 

|/|=fc m£l 

n0 



1 



2s S (v m - V n ) 



1/2 



I I exp{Tia^$d Vm ) ] [ 



me/ 



me/ 



2s S (v m - V n ) 



1/2 



(2.145) 



(2.146) 



Then the coefficients are positive and real-analytic on G. 
Next, we define a dual Toda weight function 

W(v) = | [ G(vj - v k + ia)G(-Vj + v k + ia). 

l<j<k<N 

Using the reflection and difference equations for the G-function, we obtain 

W(v) = I [ 4& + (vj - v k )s-(vj - v k ). 

l<j<k<N 



(2.147) 



(2.148) 



Thus W is entire in v and positive on G. Taking positive square roots, it is now not hard 
to verify that on G we have 



AM = W (v)~ 1/2 Hi >s (v)W (v) 1/2 , ±1 



,N, 5 = +,-. 



(2.149) 



Indeed, we have chosen the phases in f)2.132p such that these relations hold true. Note that 
the AAOs Hi t s(v) f )2.146p are formally positive operators on L 2 (G,dv), so that the same 
is true for the operators Ai } g(v) on L 2 (G, W{v )dv ). Furthermore, we have the following 
obvious corollary of Theorem 12.101 
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Corollary 2.11. For any I G {±1, . . . , ±iV}, 5 G {+, — } and a G {±1}, we have the 
dual kernel function identities 

(H l>s (v) - H^s(w))W(v) l/2 W(w) 1/2 S(v, wY = 0. (2.150) 

Just as for the hyperbolic case, we have found kernel functions relating dual Toda 
AAOs Ai t s(v) in N variables vi, . . . , to dual Toda AAOs in N— £ variables wi, . . . , %_< 
for any £ £ {1, . . . , N — 1}. However, in this case there appear to be no kernel functions 
for both signs of i at once. Another difference with the hyperbolic case is that the kernel 
identities in the following theorem and its corollary only involve two AAOs. 

Theorem 2.12. Define kernel functions 

(ri £ ( N N ~ £ W N N ~ £ 

2a - ( _ ^2 w n ) ) II II G ( Vm ~ w ™- ia )' (2.151) 

\m=l n=l J J m=l n=l 



where t = +,—, and £ = 0, 1, . . . , N — 1. For any k G {1, . . . , iV — £} and 5 G {+, — }, we 
have 

A k)S (vi, v N )SJ(v, w) = A_fc )( y(u>i, . . . , w N -i)SJ{v, w). (2.152) 
Proof. Starting from the £ = case (12.1401) . induction on £ readily yields 

r^ri An 5 8 t( v n ~ v m) eS ^ Vm ^ II s (y _ w ) 

lc{i,...,N\ mel mel mg/ 

\I\=k nil ne{l,...,N-£} 

E II ;I[e 5 (£w m ) J] — (2.153) 

|/|=fc ng/ ne{l,...,AT} 

Indeed, assuming (12.1531) for some £ G {0, 1, . . . , N — 2}, the case £+1 follows upon taking 
Wn_i to oo. 

The crux is now that the difference equations (IA.24j) for the hyperbolic gamma function 
imply 

• ■ N N-l N-l 1 
% \ -n- 1 



x s s (v m -w n - ia_ S ) 



( i\ T-r 

Sl{v,w)~ l ey^{-ia_ 5 d Vm )SJ(v,w) = e s {rt(v m -ia_ 5 /^)) ( - J JJ 

(2.154) 

SJ^w^exp^a-AJSJ (?;,«;) = e 5 (r£(w m + za_ <5 /2)) (-J JJ 



• v JV jv .. 

* \ TT 1 



(2.155) 

Hence, after shifting the variables t> m to v m + ias/2 and the variables w n to w n — ias/2, 
it is readily seen that the kernel identity (12. 1 52j) is equivalent to the functional identity 
(12.1531) for the case r = +. Taking v, w — > —v, —w in (12.1531) . we can pull out the signs 
from the s^'s so as to obtain the functional identity equivalent to the kernel identity 
(I2TT52D with r = -. □ 
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As a corollary, we shall now obtain two more kernel functions Sf(v,w). First, since 
( 12. 152ft only involves AAOs of the same order k, the identities remain valid if we multiply 
Sg(v,w) by a factor 

2 N 



N N-t 



exp [Y] v m ~ V w n \ \. (2.156) 

^ 2a + a_ ^ fx J J 

By using (I2.135[) we thus obtain 

A- k)S (vi, v N )SJ(v, w) = A k ,g(wi, w N -e)<SJ(v, w), k = 1, . . . , N — £, (2.157) 
where 

N N-t 

-I' 



SJ(v,w) = expl- Tm [^2v m -^2w n \ \SJ(-v,-w). (2.158) 



\m=l n=l 




Moreover, it is readily verified that the additional exponential factor in these kernel func- 
tions entails the identity 

A N>s (v u . . . , v N )SJ{v, w) = A-( N -g) iS (wi, • • • , w N _ e )SJ(v, w). (2.159) 

We now act with the A AO A N>s (v) on (l2~T57j) . and then use (l2~T59l) and (12TT37D . Finally, 
taking — k — > k, we obtain the following corollary of Theorem 12.121 

Corollary 2.13. We have eigenfunction identities 

Ae jS (vi, v N )SJ(v, w) = SJ(v, w), (2.160) 

and kernel identities 

A kjS (vi, . . .,v N )SJ(v,w) = A_( fc _^) j5 (wi, . . . ,w N -e)SJ(v, w), k = £ + 1, . . . , N. (2.161) 

We note that the four kernel functions Sf (v , w) and Sf{v, w) are closely related to the 
function Q in [KLS02j . which plays an important role in a construction of eigenfunctions 
of the relativistic periodic and nonperiodic Toda systems, cf. Eqs. (3.46), (2.8) and (A. 23) 
in |KLS02| . 



3 Backlund transformations 
3.1 The elliptic case 

As explained in Section 1, we are going to study the 'expected' classical asymptotics ( 11.311) 
of the elliptic kernel function ^(x,y) (\2.8\\ by introducing dependence on h in a quite 
special way. It is in fact easy to see how this should be done, because we have started 
our account in Section 1 with a description of the quantization of the classical systems 
that involves h explicitly. Specifically, we need only revert from our parametrization of 
the two positive step sizes a + and a_ in the elliptic gamma function to the parameters a 
and hp, cf. fll.23p . Taking h to then amounts to taking a_ to 0. (We keep f3 = 1/mc 
fixed, since we wish to stay in the relativistic setting.) Therefore we can study (11.311) via 
the limit (1A.14|) . Recalling the definitions of the functions W(z) and S(x,y) featuring in 



^(x,y) (cf. (!2T9|) - (l2~T2|) ). the following lemma easily follows from ( ]A.4j) and (1A.14I) . 
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Lemma 3.1. Let p G i(0,a). For x and y in the elliptic configuration space G (jl.8p . we 
have classical limits 

\imihlnS(r,a,h(3]x,y) = — 2_] / diuln R(r, a; w), (3.1) 

I r'Xj-Xk+ia/2 

limiMn W(r, a, h/3;x) — — } / dw In R(r, a; w), (3.2) 

^° P Jxi—Xh+ia/2—p 



where the integration paths stay away from the cuts given by f ]A.13j) . 

On account of the restrictions x,y G G and p G i(0, a), we can actually use the 
representation ( ]A.15j) for lni? on the right-hand side of (13. ip and ( 13. 2p . 

To begin with, we now analyze whether the function F(x, y) resulting from the above 
limits according to (ll.3ip gives rise to a Backlund transformation. Thus, we study a trans- 
formation B from the canonical coordinates (x,p) G Q (with the elliptic phase space Q 
given by (11.7p -( fL~8"]l ) to new canonical coordinates (y, q), by taking as generating function 

F(x,y) = ^(F w (x) + F w (y)+F s (x,y)), (3.3) 

where 

Y rXj-x k +ia/2 

F w(x) = -J2 dw\nR(r,a;w), (3.4) 

J x 3 -x k +ia/2-p 

N rxj-yk~p/2 

Fs(x,y) = / dw In R(r, a; w). (3.5) 

By definition, this means that y(x,p) is to be determined from the equations 

dF 



Pj 



dxj 



1 j f R(xj — x k — ia/2)R(xj — Xk — p + ia/2) 



2 P{^ \R(x j -x k + ia/2)R(x j -x k + p-ia/2)J (3.6) 

N 

/ 1(1 .;• : — //,. 1)1 v s \ 

j = l,...,N, 



If, ( R( Xj -y k + p/2) \ 



and then q(x,p) is given by 
dF 



-Vln 
7R L^i 



R{Vj ~Vk + ia/2)R(y j - y k + p- ia/2) 



~ 2(5 f£ \R( yj -y k - ia/2)R{y 3 - y k - p + ia/2) J (3.7) 
(We used evenness of lni?(z) in these formulas, cf. (IA.15p .) 
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We have now arrived at the point where we can elaborate on the problems alluded to 
in the Introduction. Ideally, we would like B to define a bijection on the elliptic phase 
space Q. Now whenever an arbitrary real- valued function F(x, y) is used as a generating 
function for a map on the phase space Q, local canonicity is clear (at least, when one takes 
for granted that there exists a local solution y to the implicit equations (13. 6p ). but in 
general it will not yield a global symplectomorphism of Q. For a special choice of F(x, y), 
therefore, it may be intractable to prove that it does. In the case at hand, however, this 
question is easily decided negatively when one retains the parameter restrictions we have 
imposed: Assuming that for a given (x,p) G Q there exists a solution y to the system 
of equations (I3.6p . this solution cannot belong to the elliptic configuration space G (II. 8p . 
Indeed, assuming y e G, it is immediate from (13 .6p that the numbers Pi, ■ ■ ■ ,Pn are purely 
imaginary, a contradiction. 

To try and get around this snag, it appears inevitable to require that (3 be purely 
imaginary instead of positive. (The requirement that p be real instead of purely imag- 
inary still leads to momenta that are not real; cf. also (ll.ll) - (ll.3l) to see why a real p 
is troublesome.) Before analyzing this change in some detail, it is expedient to study 
first in what sense the map might be a Backlund transformation. Reasoning formally 
(in particular, assuming its existence at least for unspecified parameters and phase space 
variables), this is readily answered by using (I3.6P and (13.71) : These equations do imply 
the Backlund property 

S k (x,p) = S k (y,q), k = l,...,N, (3.8) 

where the Hamiltonians are given by (ll.ll) - (ll.3p . We proceed to explain this. 

First, we note that in view of flA.20p and (IA.12[) the interaction function (11.31) can also 
be written 

/(zj-exp( irp)^ R{z _ ia/2)R{z + ia/2) J • (3.9) 
Hence (13. 8p is equivalent to 

\ 1 /2 

ETT i a \ TT / R ( x m ~ %n + P - ia/2)R(x m - x n - p + ia/2) x "~ 
[[exp{(3p m ) [[ 

IC{1,...,N} m£l mel 
\I\=k n t l 



I C{1, N} m€l m€l 
\I\=k 



R(x m — x n — ia/2)R(x m — x n + ia/2) 
R(y m -Vn + P- ia/2)R(y m - y n - p + ia/2) \ 1/2 



R{Vm -Vn- ia/2)R(y m -y n + ia/2) 



(3.10) 



Next, from (13. 6p and (13. 7\} we have 

exp((3p ) - TT ( R ^ Xm ~ Xn ~ ia / 2 )R( x m -x n - p + ia/2) \ 1/2 -pr R(x m - y n + p/2) 
~i \R{x m -x n + ia/2)R(x m -x n + p-ia/2)J ^ R(x m - y n - p/2)' 

(3.H) 

(a \ — TT ( R (y™ ~Vn + ia/2)R(y m - y n + p - ia/2) \ 1/2 -p-r R(x n -y m + p/2) 
ex P {pq m ) - ll _ ^ _ ia/2)R{ym _ yn _ p + ia/2) J 11 R(Xn _ ym _ p/2) ■ 

(3.12) 
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Consider now the product of the quantities exp((3p m ) for m in a fixed index set /. For 
a pair mi, 777.2 G ^> the two corresponding radicand terms coming from the first product 
in (13. lip cancel, since R(z) is even. For pairs m G I, n £ I, we can cancel two of the four 
radicand factors and combine the remaining two to rewrite the left-hand side of (I3.10p . 
Likewise, the right-hand side of (I3.10P can be simplified. The upshot is that the Backhand 
property (13. 8p holds, provided the following identities are valid for k = 1, . . . , N: 

y- t-t R(x m -x n - p + ia/2) yr R(x m -y n + p/2) 

| 7 | =fc n0 n=l,...,N 

\- t-t R(y m -y n + P -ia/2) -j-r R(x n -y m + p/2) 

~ ^hmlh R(ym-yn-ia/2) 11 R(x n -y m - P /2y { - 

\I\=k n$I n=l,...,N 



These functional equations can be reduced to (I2.19P by using (1A.20|) ; cf. also Eqs. (2.8)- 
(2.9) in [Rui06j. Therefore, we have now demonstrated (13. 8ft . 

Next, recall that the kernel property is not spoiled when all coordinates x\,. . . ,xn 
are translated by £ and when ty(x, y) is multiplied by a function of the form (ll.29p . In 
particular, we can allow a multiplier 

N 

exp(i^(s(x,y))/h) } s(x, y) = ^(xj - yj), (3.14) 

3=1 

where ip(z) is any entire function. Taking the classical limit as before, we obtain a 
generating function that yields an extra term ip'(s(x, y)) on the right-hand sides of (13.61) - 
( 13. 7p . Since these extra terms are equal and do not depend on the index j, they do not 
spoil the argument leading to the validity of ( 13. 81) . Hence a large family of Backhand 
transformations arises. 

Returning to the non-rigorous status of these developments, we first repeat that we 
must switch to a parameter /3 that is purely imaginary to ensure that for real positions 
and momenta x,p there might exist vectors y, q that not only satisfy ( I3.6p -( l3~7|) . but are 
also real. As they stand, the Hamiltonians Sk{x,p) (II. ip are then not real- valued on Q, 
but this is easily remedied by switching to 

H k (x,p) = S k (x,p) + S k (x, -p), k = l,...,N. (3.15) 

Now this looks satisfactory at face value, but in fact problems remain. The point is that 
it seems very unlikely that the commuting local flows generated by Hi, ... , H^-i extend 
to global flows (the ifjv-flow is of course global). For H\, for example, the dependence on 
Pj is now via a factor cos(\ /3\pj). Hence, conservation of Hi is compatible with collisions 
after a finite time, and constant-ifi hypersurfaces are not compact, in contrast to the 
positive-/? case. 

Put differently, one should not expect to obtain Liouville tori in Q for (3 purely imag- 
inary. Moreover, a proof of existence and uniqueness of a solution y in the elliptic config- 
uration space G (11.81) to the equations (13.61) for a given (x,p) G Q is not in sight. Finally, 
one has to deal with the Poisson commuting Hamiltonians (I3.15P on Q, which can yield 
singularities in finite time, hence an unphysical behavior. 
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Ignoring these problems, we continue by tying in the above generating functions with 
the ones that can be found in [NRK 96J . To this end we first view the generating function 
F(x,y) given by (I3.3[) ( T375T) as the Lagrangian of a 'discrete-time' map. (In a general 
setting, the relation between these two viewpoints has been clarified by Veselov |Ves91j .) 
For a given sequence of vectors z{n) G C^, n G Z, and a fixed value of n, we write 
z = z(n), z = z(n + 1) and z = z(n — 1), in accord with the notation used in [NRK96J. 
To make the connection to [NRK96J, we shall actually start from a slightly more general 
generating function, namely, 

F t,v( x > V) = F ( x i +€,---,x N + £,y 1 ,...,y N ) + 7s(x, y) 2 //3, (3.16) 
cf. the paragraph containing (|3.14p . Then from the discrete Euler- Lagrange equations 

^M + ^M) = 0, j = l,...,N, (3.17) 

and (I3.6p -( l3~7|) we deduce 

R(zj - z m — p + ia/2)R(zj — z m — ia/2) 
, . R(zj - z m + p - ia/2)R(zj — z m + ia/2) 



Pin) II 



R(z m - Zj + p/2 + i)R{z 3 - z m - p/2 + 
= M R(z m - Zj ~ p/2 + 0R( Zj -z m + p/2 + ' (3 ' 18) 

where the prefactor reads 

JV 

P( 7 ) = J] exp(2 7 U m + z m - 2z m )). (3.19) 

m=l 

Next, we choose 

C = p/2 + ia/2, (3.20) 
and use (1A. 12|) and (1A.20P to rewrite the result as 

P( 7 ) TT s J^i ~ Z ™~P) = TT s{zj-z m -p)s{z j -z m ) ^ = ^ ^ ^ 

^j s ( z j- z ™ + p) - z m )s(z j - z m + P y 

We are now in the position to compare f)3.2ip to Eq. (2.14) in |NRK96] . To this end 
we substitute 

p = -A, (3.22) 

and use (IA. 18[) to switch from the s-function to the Weierstrass cx-function. Then (13.211) 

becomes 

pd - * a) n :[r t + » = n :ir M t 1 1 x 

<J\Zj — Z m — A (7{Zj — Z m CT{Zj — Z m — A 

mjtj J m=l J J ' 

(3.23) 

This coincides with Eq. (2.14) in |NRK96j . provided p/p is given by Eq. (6.20) with 8 = 0, 
and 7 is suitably specialized. 
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Even though (13. 2 1 p can be made to coincide with the 'discrete-time Newton equations' 
(2.14) in |NRK96] . the generating functions and Lagrangians employed in |NRK96] differ 
from the above ones. This is due to a different choice of phase space variables implicit 
in |NRK96] . which leads to their functions having an asymmetric dependence on x and 
y. (The choice amounts to a canonical map of the form (x,p) (x,p+ f(x)).) 

We should add that the existence of sequences of vectors z(n) G satisfying (I3.2ip 
is left open. In any case, it seems unlikely that for given initial values z(0),z(l) in the 
elliptic configuration space G (11.81) there exists a solution sequence that stays in G. 

3.2 The hyperbolic case and its dual 

From now on we reparametrize the two length scales a± as 

a+ = 2vr//i, a_ = h/3, (3.24) 

except in the elliptic cases considered in Section 4, where we retain an imaginary period 
ia. (Recall we motivated this change at the end of the Introduction, cf. the paragraph 
containing (ll.33p .) It is also convenient to trade the parameter p (which has dimension 
[position]) for a (dimensionless) parameter 

T = -inp/2. (3.25) 

With these changes, we obtain from ( 1A.29I) the following counterpart of Lemma 3.1. 

Lemma 3.2. Let r G (0, 7r). For x and y in the hyperbolic configuration space 

G hyp = {x G R N | x N < ■■■ < x-l}, (3.26) 

we have classical limits 

I N rK x j-Vk)-ir 



\imih\nS(2ir/n,Hp;x,y) = — > / dw ln(2 cosh(w/2)), (3.27) 

limihlnW{2n/p,h(3;x) = —J2 ln(2 cosh(w/2)), (3.28) 

m PL 1 TltJu(x,-X k )+i7T-2iT 



j k=1 J /j,(xj-y k )+iT 

I Pfl(Xj-X k )+ilT 

' fi(xj —x k )+iir—2iT 

where the integration paths stay away from the cuts ±i[ir, oo). 
With the change Q — > i\ yp , where 

fihyp = {(x,p) G M. 2N I x G G hyp }, (3.29) 



the developments leading to (13.151) now apply with straightforward adaptations. In par- 



ticular, the asymptotics (ll.3ip entails generating functions 



F{x,y) = ^-(F w {x) + F w {y) + F s {x,y)), (3.30) 



where 

I PH{Xj-X k )+iTT 



1 Pf-^K^j 

F w (x) = -J2 dtuIn(2cosh(u;/2)), (3.31) 

^ j-J.fr J V(Xj-X k )+iTT-2iT 
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F s (x,y)=J2 dwln(2cosh(w/2)), (3.32) 

„■ j„_i J fi{xj-y k )+iT 



j,k=l ■ 

and the Backlund property (I3.8P is equivalent to the identities 

tt sinh(/i(a; m - x n )/2 - it) y-r cosh((/i(x m - y n ) + it)/ 2) 
ic{h.,N}mei sin Hv( x m-x n )/2) 11 cos h((n(x m - y n ) - ir)/2) 

\I\=k n n=l,...,N 



ic{b.,N}ma ^HKVm-Vn)/2) ^1 cosh((/x(a; n - y m ) - it)/2) ' 

|j|=jfc n^/ n=l,...,JV 

Next, starting again from a modified generating function (I3.16p . we get from the 
discrete Euler-Lagrange equations (13. 1 7[) as the counterpart of (I3.18P 



^(7) n 



smh(fj,(zj — z m )/2 — it) 
/;j smh(n( Zj -z m )/2 + it) 



JL cosh(/i(z m - Zj + 0/2 + it/2) cosh(fi( Zj - z m + 0/ 2 - ir / 2 ) 

1J- cosh(/x(£ m - Zj + 0/2 - ir/ 2) cosh(//(zj - 5 m + 0/2 + «V/2) ^ ^ 

Choosing 

£ = zr//i + in/fi, (3.35) 

this becomes 

p il) TT sinh ^ /i ^ j ~ ' 2m )/ 2 ~ ir ^ - TT sinil (l J '( z 3 ~ &m)/ 2 ~ ir ) smh(fj,(zj - z m )/2) 
sinh(yu(zj — z m )/2 + it) "-jf- sinh(/^(£j — z m )/2) sinh(//(^- — z m )/2 + it) 

(3.36) 

This set of equations amounts to the hyperbolic specialization of Eq. (2.14) in [NRK96J. 

Just as in the elliptic case, it is unlikely that any solution sequences z(n) exist that 
stay in Gh yP for all 'discrete times' n G Z. To be sure, the equations to be solved do not 
involve 0, so it might appear that the issue whether /3 is real or imaginary is moot. In 
fact, however, the difference is decisive, since only in the latter case the equations have a 
chance to correspond to a canonical map on £\ y p- But since the commuting flows are not 
global for /3 imaginary, there exists no well-defined action-angle map, by contrast to the 
case G (0, oo) [Rui88j . 

In fact, the existence of solution sequences in C N has not been shown beyond doubt 
even in the hyperbolic case. In |NRK96] there are no reality conditions specified, and 
although the arguments for the hyperbolic case are formally convincing, they involve 
tacit assumptions that are not checked. 

Finally, we discuss the dual hyperbolic case. This can be quite easily handled after 
our notation change fl3.24l) - (l3.25l) . Indeed, we need only substitute 

x^p, p^x, y->q, q-+y, (3.37) 

replace Ghy P by 

G hyp = {pER N | pjv< •••<&}, (3-38) 
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and interchange 

<B- fi, (3.39) 

wherever these parameters occur |Rui88] . (See also Section 4, where we have occasion 
to say more about the self-duality of the relativistic hyperbolic case.) In particular, note 
that the resulting generating function F(p, q) now gives rise to purely imaginary positions 
x and y unless /i is required to be purely imaginary. 

3.3 The periodic Toda case 

We proceed to obtain the classical asymptotics ( 11.311) of the periodic Toda kernel functions 

^(x, y) = U(x) Tl/2 U(y) Tl/2 S ± (x, y). (3.40) 

We recall that U(x) and S ± (x, y) are given in terms of the functions Gr and Gl, cf. (12.711) 
and (I2.87l) -( l2~88j) . and that we have switched to the parameters (13. 24ft . Therefore, taking 
h to amounts to taking a_ to 0, so we infer from ( 1A.29I) and ( 1A.31I) -( I"A.32I) that we have 



1 f z 

\imih\nGii(2Ti/ fi, hp 1 ; z) = — / dw In (l + exp(=)=/iw)) . (3-41) 

R4-0 L p Jo 



The following lemma is now easy to verify. 
Lemma 3.3. For x and y in M, N we have classical limits 

rx m+1 -y m ±i-K/2n+ri/2 



limi/Un S ± (2n/ /i, h(3; x, y) = \ J dw In (l + exp(fiw)) 

/ym-Xm±iTr/2ii+ri/2 \ 
dw In (1 + exp(fiw)) J , (3.42) 



\imih~\n U(2tt/h, hp 1 ; x) = — — / dw In (l + exp(/iu>)), (3.43) 

m=l 

where the integration paths stay away from the cuts ±z[?r//i, oo). 

Taking fl 1 . 3 1 [) as a guide, we should consider transformations B ± from canonical co- 
ordinates (x,p) G M. 2N to new canonical coordinates (y,q), generated by 

FHx,y) = ^-(F u (x) + F u (y)-F^(x,y)), (3.44) 



where 



1 N rx m+1 -x m +r) 

Fu(x) = -y2 dw1n(l+exp(jiw)), (3.45) 

2 m=X Jo 



(rx rn+ -i-y m ±iTT/2fi+ri/2 
/ dw\n (l + exp(/iw)) 

Jo 

i>ym-Xm±iTr/2fi+r)/2 \ 

+ J dwln (l + exp(/iw)) J . (3.46) 
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However, as announced in the Introduction, there is a problem with these generating 
functions in the Toda well. Indeed, we need r\ to be real to ensure formal self- 

adjointness of the Toda AAOs, but this requirement implies that the gradients of F + 
and F~ are not real- valued on ~R 2N . Hence we should not expect that they give rise to 
bijections on the Toda phase space M? N . 

To attempt to improve this state of affairs, let us first focus our attention on F + (x, y) 
and consider a suitable analytic continuation. Specifically, taking 

x m -> x m - irc/2fi + r]/2, m = l,...,N, (3.47) 

the function Fjj(x) is left invariant, while F^(x,y) turns into 



N 



»Xm+l—ym+ri 



m=l 



F s( x > V) = ( / dtw In (1 + exp(fiw)) 

dw In (l - exp(nw)) J + NI/fi, (3.48) 



where / is the integral 



dzln(l + exp(z)). (3.49) 



Since additive constants are irrelevant, we have now been led to a more well-behaved 
generating function. Clearly, if we start from F~(x,y) rather than F + (x,y), we would 
obtain the same function up to an overall sign (which is a matter of convention anyway). 
The key question remains, however, whether this function does give rise to a symplecto- 
morphism (x,p) i— > (y,q) on M. 2N via f ll.30p . Accordingly, reverting from now on to the 
dimensionless coupling 7 defined by (recall (11 . 16[) and (11 . 22[) ) 



7 = exp(/x///2), (3.50) 



we should study the equations 

dF+ 



Pr, 




1 + 7 2 exp (fi(x m+ i - x m )) 

1 + 7 2 exp (n(x m - x m _i)) ) ( 3 ' 51 ' 
1 + 7 2 exp (fx(x m - y m -i)) 



1 - exp(n(y m - x m )) 



771=1 N, 



dF+ 




1 + 7 2 exp (fj.{y m - y m -i)) 

l + 7 2 exp (//(y m+1 - y m )) ) ( 3 ' 52 ) 

1 + 7 2 exp (fj,(x m+ i -y m )) . 
1 - exp (^(y m - x m )) 
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There is still a troubling aspect about these equations: If we start from an arbitrary 
phase space point (x,p) G R 27V , then we can only accept a solution y(x,p) to the implicit 
equations ( 13. 51 f) when it satisfies 

y m < x m , m = l,...,N. (3.53) 

Now the existence of a solution with this property cannot easily be ruled out, but also 
seems hard to prove. Even when we assume that such solutions exist and give rise to a well- 
defined symplectomorphism, there is still a problem with reinterpreting the generating 
function as a Lagrangian for a discrete map on ~R N . Indeed, in that case we would expect 
to have the freedom to choose arbitrary x, y in M. N to obtain a unique sequence z(n) G R , 
n G Z, with initial values z(0) = x, z(l) = y, whereas this freedom is at variance with the 
constraint (I3.53p . 

Ignoring these dilemmas from now on, we can easily show the Backlund property (13.81) 
by proceeding in the same way as in Subsection 3.1. Specifically, in this case (13. 81) can be 
written 



+ 7 2 exp (n(x m+ i - x m )) ) 
Y[ (l + 7 2 exp (/i(x m -x m _i))j 1 = (x,p->y,q). (3.54) 



IC{l,...,N}m£l m€l 
\I\=k m+lil 



Also, from (l33Tj) -( l3~52]) we have 

1 /2 

, Q s ( l +7 2 exp (n(x m+1 - x m ))\ l + 7 2 exp(/i(x m -y m _i)) 

exp((3p m ) = —— — — — ^— , (3.55) 

y 1 + 7 2 exp [fjt{x m - Xm-!) ) J 1 - exp [fj,{y m - x m ) ) 

exp m q ) - ( 1 + 7 2 exp (fj,(y m - y m -i)) \ 11 ~ 1 + 7 2 exp (fi(x m +i - y m )) ^ ^ 
y 1 + 7 2 exp (n(y m +i - y m )) j 1 - exp (fi(y m - x m )) 

Now fix an index set I and consider the product of quantities exp(/3p m ) with m in /. 
For m\,m2 G / such that m<i = ni\ + 1 (mod N), the denominator of the radicand 
corresponding to m,2 cancels the numerator of the radicand corresponding to mi, and vice 
versa for m2 = mi — 1 (mod N). The product of quantities exp(/3q m ) can be simplified 
similarly. From this it readily follows that the Backlund property f l3.54j) amounts to the 
functional equations 



1 + 7 2 exp (y.{x m - y m -i)) 



hit Xuijj 



nf, , 2 ( ( \\\ TT i + 7 exp [fi[x, 

(1 + 7 exp (/i(x m+1 - x m )) j || - 

+10 

\T TT /, . 2 / / . A TT ! + 7 2 exp (/i(x m+ i - y m )) 

-. meI meI 1 exp {(j,\y m x m )j 

\I\=k m-HI 



\I\=k rn+10 



Ic{l,—,N] m€l m€l 



These equations can be reduced to the identities ( 12.571) by taking 

5 = +, a + = 27r//x, p = (2 In 7 + in)/ fi. (3.58) 



41 



Thus we have established ( I3.8P in the periodic Toda case. 

Finally, we tie in the above generating function with a 'discrete-time' map introduced 
by Suris |Sur96 ] as a time-discretization of the defining Hamiltonian of the relativistic 
periodic Toda system. To this end we observe that the discrete Euler-Lagrange equations 

— (z,z) + —(z,S) = 0, m = l,...,N, (3.59) 
oy m ox m 

for a sequence of vectors z(n) £ C N , n £ Z, are equivalent to 

1 - exp (fj.(z m - z m )) _ 1 + 7 2 exp (fi(z m+1 - z m )) 
1 - exp (n(zm - z m )) 1 + 7 2 exp (fx(z m+1 - Z m )) 

l+7 2 exp(/i(2; m -5 m _ 1 )) 

x — — )— rf, m = l,...,iV. (3.60) 

1 + 7 exp (/i(2; m - z m _i) J 

These equations coincide with Eq. (5.4) in |Sur96] if we substitute fj, — 1,7 = g. 
3.4 The nonperiodic Toda case and its dual 

The classical asymptotics (jl.3ip of kernel functions (13.401) with U(x) given by (12.991) and 
S ± (x,y) by (12. 101|) (12. 102j) can be studied in virtually the same way as in the periodic 
case. The reader who has followed us this far will have no difficulty to make the pertinent 
changes. In particular, with the nonperiodic Toda convention (11. 18ft in force, the Backhand 
property ( 13.81) continues to hold, and the nonperiodic analog of F + (x, y) can be readily tied 
in with Suris' time-discretization of the defining Hamiltonian of the relativistic nonperiodic 
Toda system [Sur96] . For this case the action-angle map is known from [ Rui9 0j , and so 
a further study might clarify the puzzling constraints (13.531) . which arise in the same way 
as in the periodic case. This is beyond our present scope, however. 

Continuing with the dual system, we get from Subsection 2.5 a kernel function 

*(p, q) = W^p/^W^q/^S^p/^ Pq/»), (3.61) 

where W and S are given by (12.1481) and (12.1331) . Also, as explained earlier, we use the 
parameters ( 13.241) . The classical asymptotics ( 11.311) is now easily inferred from ( 1A.29I) . 



Lemma 3.4. For p and q in G ( 12.1451) . we have classical limits 

N r /3(p m -q n ) 



\imihlnS(27r/fJ,,h/3;Pp/fji,/3q/fj,) = — V / dw In (2 cosh(w/2)), (3.62) 

r4 ° m,n=l J ° 



]xmihhxW(2Tr/fi, h/3; fip/fi) = — ^ (Pm-Pn), (3.63) 

r 1< 



'171 

l<m<n<N 

where the integration paths stay away from the cuts ±i[7T, oo 
Hence we should study whether the generating function 



F(p,q) = ^{F w (p) + F w (q) + F s (p,q)), (3.64) 
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where 

N 
m=l 

N fP{pm-qn) 

F s {p,q)=22 / dwln (2cosh(w/2)), (3.66) 

i Jo 

m,n=l " 

might lead to a canonical map (x,p) (-> (y, q) on the dual Toda phase space 

Q = {(x,p) eR 2N \peG}. (3.67) 
This entails that q(x,p) is to be determined from the equations 

i 1 N 
x m = — 2m + 1) In (2 cosh(/3(p m - g„)/2)), (3.68) 

and then y(x,p) is given by 

iV 



27T 

Vm 



(jV-2m + l)--^ln(2cosh(/3(p n -g m )/2)). (3.69) 



71=1 



Clearly, for p,q £ G the numbers x m and y m are not real as they stand, and this 
cannot be cured by switching to a parameter p, that is purely imaginary. On the other 
hand, we can also start from a modified generating function 

JV 

F(p,q) = — ^2(p m -q m ) + F((p 1 -m//3 } ...,p N ~iTr//3) } q) } (3.70) 
Zfl m=l 

which can be obtained from a modified kernel function, as explained earlier (cf. the para- 
graph containing (jl.29p ). The point of this is that a solution q G G to the modified 
equations 

dF m 



y (iV-2m + 2)-i^ln(-2zsinh(/3(p m -g n )/2)), (3.71) 
Op m Zfi p 

might then exist, provided its coordinates interlace with those of the given vector p G G: 

q N <p N < q N ^ < ■ ■ ■ < qi < p 1 . (3.72) 
In that case, the vector y given by the modified equations 

dF 1 N 

y m = ^~ = ^-(N - 2m) - - Vln ( - 2% sinh(/3(p n - q m )/2)), (3.73) 
dq m 2/i ^t^i 

would be real. 

Again, a further study of this global analysis question is beyond our scope. Rather, 
we proceed to show the validity of the Backhand property 

H k (x,p) = H k (y,q), k = l,...,N, (3.74) 
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for the dual Hamiltonians, which we rewrite as 

P/2 



H k (x } p) = Yl Uv~ 

t J. J. Oil 



. sinh(/3(» m — p n )/2)\ 

/C{1, .. .,iV} mel 1 V^V^m mil j\ leJ 

\I\=k n$I 



JTexp(//xj), k = l,...,N. (3.75) 



Substituting the above expressions fl3.68p - fl3.69p for x and y, we infer that f)3.74p holds 
if and only if the following identities are valid for k — 1, . . . , N: 



IC{1,...,N} m£l 

\I\=k 



e ri ex p(-T( iv - 2m+i ))n 

n<£I 



1 



n 



sinh(/3fj3 m -p n )/2) 
1 



m£l 
rc=l,...JV 



cosh(/3(p m - q n )/2) 



IC{1,...,N} m£l 

\I\=k 



e n ex p(T^- 2m+i ))n 



n 



sinh(/3(g m - g n )/2) 
1 



n=l JV 



cosh(/3(p n - q m )/2) 



(3.76) 



Fix an index set / and consider the corresponding term on the left-hand side. Since p 
belongs to G f)2.145p . we can rewrite this term as 



n 



-JV+2m-l 



iic-mi 



N 

n 



i 



m£l 



nil nil Shl HP(Pn ~ Pm)/2) COsll(/3(j3 m - g n )/2) ' 

n>m 

The numerical factor can be simplified using 



(3.77) 



n<->= n H/ n<-)=ni 

m£l m£l I rn,n£l m£l 

n£I n=l,...,N n>m 

n>m n>m 



viV-m / /_\fc(ft-l)/2 _ A~k(k 



1 ii 1 



2N-2m 



Doing so, we obtain 



JV 

n 



11 sinh(/3(p n -p m )/2) IJ^ cosh(/3(j3 m - q n )/2) 



(3.78) 



(3.79) 



Likewise, since q G G, the corresponding term on the right-hand side can be simplified to 
yield 

1 N 1 

MN-k) TT L _ TT L _ ,o on) 

11 sinh(/3(g m - q n )/2) 11 cosh(/3(p n - q m )/2) ' 1 ' ) 



We can thus reduce (I3.76P to (I2.140p by substituting pj — > pj — nr//3, j = 1, . . . , N, and 
cancelling the /-independent numerical factor. This proves the Backlund property f)3.74p . 
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4 Nonrelativistic limits 



4.1 Kernel functions 
4.1.1 The elliptic case 

To take the nonrelativistic limit c — > oo of the various quantities in Subsection 2.1, we 
first need to reparametrize the two positive step sizes a + and a_ in the elliptic gamma 
function as a and h(3, cf. ( 11.231) . Taking /3 = 1/mc to then amounts to taking a_ to 
0. (We keep h fixed, since we wish to stay in the quantum setting.) To obtain nontrivial 
limits, we should also set 

p = i{3g. (4.1) 

Since we let go to 0, the coupling constant g is allowed to vary over (0, oo). 

From ( 1A.16I) it now follows that the kernel function S(x,y) (12.12ft has — > limit 



N 

K{x,y) = H R{x 3 -y k )- 9 ' h . (4.2) 
j,k=i 

Moreover, ( 1A.16I) yields the nonrelativistic limit of the weight function given by ( 12.91) - 

] [ R(xj - x k + ia/2)R(xj - x k - ia/2) J . (4.3) 

l<j<k<N J 

The nonrelativistic counterparts H miQr , m = 1, . . . , N, of the commuting Hamiltonians 
H k .+ (12.11) arise as — > limits of suitable linear combinations of the Hamiltonians 
. . . , H m ^ + . This limit transition hinges on the use of the classical elliptic relativistic 
and nonrelativistic Lax matrices. Later on, we use these matrices in our study of classical 
Backhand transformations, cf. (14.1311) and (14.1331) . On the quantum level, however, the 
pertinent limit is fraught with ordering problems. Although these can be resolved, the 
details are quite substantial and will be skipped. They can be found in Subsection 4.3 
of [Rui94j . (We employ a similar method in the quantum periodic Toda case below, and 
make use of a corresponding formula for the Hamiltonians in the hyperbolic case discussed 
in the next section.) 

An important ingredient of the reasoning in |Rui94] is a uniqueness result obtained 
by Oshima and H. Sekiguchi [OS95] ; this paper also contains explicit expressions for the 
commuting elliptic PDOs. The relevant limits entail the kernel identities 

{H Km {x) - H k ^{-y))m m {x, y) = 0, k — 1, . . . , N, (4.4) 

where 

y ni (x,y) = W nr (x) l / 2 W nr (yf' 2 lC{x,y). (4.5) 
More explicitly, the Hamiltonians are commuting PDOs of the form 

N 

Hx, m (x) = -ihJ2 d Xi, ( 4 - 6 ) 
H 2 ,nr(x) = -h 2 ^ d x 01 d X]2 - g(g - h) ^ P( x j ~ xi]ir/2r,ia/2), (4.7) 
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H k>aT (x) = (-ih) k d Xji ---d Xjk +\.o., k = 3,...,N, (4.8) 

l<ji<-<Jfc<JV 

where 1. o. denotes terms that are of lower order in the Xj-partials [OS 95]. (In ( 14. 7ft we 
have omitted an additive constant that depends on the spectral parameter in the Lax 
matrix.) In particular, for the defining Hamiltonian 



l h 2 N 

Hnr(x) = -(Fi inr (x)) 2 - H 2 ,nr(x) = - — ^ d l 3 + 9(9 ~ ^1 P( X 3 ~ *"/ 2 ) 

j=l l<j<l<N 

(4.9) 

of the elliptic nonrelativistic Calogero-Moser system, this implies the kernel identity 

(H m (x) - H nr (y))* nr (x, y) = 0. (4.10) 

(Here and below, we choose m — 1.) This identity was first obtained by Langmann 
|Lan00] . cf. also |Rui04] for further details. 

4.1.2 The hyperbolic case and its dual 

Using the parameters (I3.24[) and (14. ip . we obtain from ( 1A.30|) the nonrelativistic limits 

N 

K(x,y) = J] pcoBh^ - y*)/2)]-"/* (4.11) 
j,k=i 

(\ sA 
II 4sinh 2 (/i(x J -x fc )/2) , (4.12) 
l<j<k<N J 

as the analogs of The kernel function ^ nr is given by (jOJ) with ([4TTT]) - (j4TT2]) 

in force. Then we obtain (I4.4[) ( ]47T0|) with the replacement 

p(xj -xi;ir/2r, ia/2) -> y u 2 /4sinh 2 (/i(x j -x{)/2). (4.13) 

In contrast to the elliptic case, we are also able to obtain kernel identities relating the 
PDOs Hk nr in N variables X — f X\ ? . . . , X TV ) to a sum of the PDOs Hj nr in N — £ variables 
y = (yi, . . . , yN-t), £ = 0, . . . , N. The key for doing so consists of explicit formulae for 
the PDOs, which involve not only the nonrelativistic Lax matrix 

L m (x,p) jk = Sjkpj + (1 - 5 jk ) . -. ifig _ — — -y, j,k = l,...,N, (4.14) 

but also the diagonal N x N matrix 

E(x) = diag(zi(x), . . .,z N (x)), (4.15) 

with 

Zj (x) = coth ~ Xfc )/ 2 ) ' 3 = h---,N. (4.16) 

We note that the substitution Xj — > ixj in L nr (x,p) yields a trigonometric Lax matrix that 
is slightly different from the Lax matrix due to Moser |Mos75] . (Specifically, in Moser's 
Lax matrix the function 1/ sin is replaced by cot.) 
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The canonical quantization substitution 



Pj ^-iHd Xj , j = l,...,N, (4.17) 

in L nT (x,p) yields an operator-valued matrix whose symmetric functions Sfc(L nr )(a;) are 
well defined, since no ordering ambiguities occur. (Indeed, a term in the expansion of a 
principal minor of (I4.14p that depends oiapj does not depend on Xj.) However, S 2 (-^nr)(^) 
is not equal to 

H 2 , nr (x) = -ft 2 d Xji d Xj2 -g(g-h) ^ //^sinh 2 ^ - x{)/2), (4.18) 

i<ii<i2<iv ~t<j<i<N 

since the term proportional to ft in the potential energy is missing. (In fact, there seems to 
be no complete proof in the literature that the N PDOs T>k{L nr )(x) commute for arbitrary 
N, cf. [Rui94j . Subsection 4.2.) 

By contrast, the symmetric functions Hk{L nr (x,p) + E(x)) with k > 1 contain products 
of terms that depend on p m and x m , so their canonical quantization is ambiguous. As 
shown in Section 4.3 of |Rui94j . the ordering choice ensuring commutativity is normal 
ordering: the procedure of putting x-dependent coefficients to the left of monomials in 
the momentum operators —ihd Xl , . . . , —iftd XN . We shall write : £fc(L nr + E)(x) : for the 
normal-ordered PDOs obtained from Efc(L nr (x,p) +E(x)) by substituting —iftd Xm for p m . 
The nonrelativistic commuting Hamiltonians Hk jQI are then given by the formula 

H k , m (x) = W^x) 1 ' 2 : % (L ni + E) (x) : W nr (x)~ 1/2 , k = l,...,N. (4.19) 

(To get a feel for what is going on here, the reader may wish to check the case N = k = 2.) 

The corresponding kernel identities have the same structure as the identities (I2.33|) in 
Theorem 12 . 11 However, in this case they involve coefficients c^j with £ e N, j G Z, given 
by 

c ,o = l, Q J= 0, j>£, j<0, (4.20) 

and 

ci,j = (^f) 3 Si(l, 3, • • • , 21 - 1), j = 0, ...,£, (4.21) 

where Sj(a\, . . . ,ag) is the jth elementary symmetric function of aj, . . . , ag. We note that 
the coefficients are uniquely determined by the recurrence relation 

ci+ij = c e>j + ^(2£ + l)c t ^ x (4.22) 

together with the side conditions (14.20)) . 

We are now ready to state and prove the pertinent identities. 

Theorem 4.1. For £ = 0, 1, . . . , N, let 



Kifa V) = TTn — nN-er n TJ-f wom /r ' N ^ = 6Xp — W 1^ 



N 



nLintipcosMM^-^)^)]^' - h 

(4.23) 
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For any k 6 {1, . . . , N}, we have 
■ S fc (L nr + E)(xx,.. .,x n ) : JC £ (x,y) 

min(fc,£) 

= ^2 Ci >j '■ ^ k -j ( Lm + E ) • • • ' ~y N ~z) '■ ( 4 - 24 ) 

i=o 

where 

:t m (L nv + E){- yi ,...,-y N _ e ) : =0, m > N — £, : ± (L m + E) : = 1. (4.25) 

Proof. We shall prove the statement by induction on £. For £ = the kernel identity 
(14.241) is equivalent to the hyperbolic counterpart of (14. 4p . Hence, we now assume (14.241) 
for £ > and establish its validity for + 1. 

Setting 

m = ^P f ff ( J] y„ - (N + <?) yjv _, ] ] , (4.26) 



n=l 



we start from the limit 



lim <f>(y u . . . , y N _i - A)JC e (x, y u . . . , yjv-^ - A) = /Q+i(x, y), (4.27) 

A— 5>oo 

which is readily verified. To make use of this, we note the commutation relation 

<f>{y):± k (L m + E){-y): 

= : % (h m + E- ^diag(l, . . . , 1, -N - £)) (-y) : <j>{y). (4.28) 



Substituting yjy-e ~ > Vn-£ — A in the Lax matrix L nr (— y), it is clear that the matrix 
elements (L nr ) N _^ k and {L nt )j^_^ where j,k = 1, . . . , N — £ — 1, vanish in the limit 
A — > oo, whereas (inr)jv-^JV-^ = ^9 yN _ t remains the same. We also note the limits 

lim z N _ e (-y u . . . , -y N _ t + A) = -^(iV - £ - 1), (4.29) 

A— i>oo 2 

lim . . . , -y N - e +A) = -z-+Zj(-y h -y^-i), j = 1, . . . , iV-^-1. (4.30) 

A— >oo 2 

We now multiply both sides of (14. 24[) by </>(?/), substitute ytq-t Vn-£ ~ A, and consider 
a term on the right-hand side corresponding to some j = 0, . . . , min(fc, £). Using the 
commutation relation fl4T2"gj) and the limits (1P7j) and (O3)) - (jO0]) . we find that the 
A — > oo limit of the term is given by 

: E k -j(L nT + E)(-yx, . . . , -y N -e-i) : /Q + i(a;,y) 

+ % -^{2£ + 1) : (L nr + £) . . . , -y N ^) : y). (4.31) 

If we now take j — > j — 1 in the sum resulting from the second term, and compare the 
result with the recurrence relation (I4.22p that uniquely determines the coefficients qj, 
then we arrive at ( I4.24[) for £ -»• £ + 1. □ 
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We continue by deducing nonrelativistic analogs of the three corollaries fl2.44p - fl2.46p . 
As a counterpart of (12.441) . we shall first consider the kernel identities involving the defin- 
ing Hamiltonian 



H nr (x) = W nr (x) l/2 



Q (: E 1 (L m + E)(x) :) 2 - : t 2 (L nr + E) (x) :J W^Or)' 1 / 2 , (4.32) 



cf. (JM]) and fQ9l . From pill we infer 



2 '■' - V 2 J \2 



1 cj 1-^,2= r-^f] (ls 1 (l,3,...,2£-l) 2 -S 2 (l,3,...,2£-r 



^P 2 (l,3,...,2£-1) (4-33) 



24 

where Pi{a\, . . . , at) denotes the second power sum symmetric function of a%, . . . , at- This 
yields the following corollary. 

Corollary 4.2. For £ = 0, . . . , N , we have 

2 2 

(H ni ( Xl , ...,x N )- H m ( yi , . . . , y N _ t )) * /)nr (x, y) = -^|-£(4£ 2 - l)^, nr (x, y), (4.34) 
where 

* t>DI (x,y) = W nr (x) 1/2 W nr (y) 1/2 JCt(x,y). (4.35) 

We note that, apart from the exponential factor in the kernel function (I4.23p . the 
kernel identity (I4.34p coincides with the hyperbolic limit of an elliptic identity obtained 
by Langmann [Lan06j . The identity can also be obtained as a special case of Corollary 2.3 
in |HL10j . More specifically, the latter corollary depends on two polynomials a and /3, 
and one function z(x), which should be fixed according to 

a(z) = z 2 , P(z) = z, z(x) = e x . (4.36) 

In addition, one should set k = g, M = N — I and N = M = 0. This results in a kernel 
identity that is equivalent to (I4.34p for H = // = 1. To be precise, the corresponding kernel 
function has a slightly different exponential factor, and the kernel identity contains an 
additional overall factor of 2. 

We proceed to detail a nonrelativistic analog of (12.451) . From (14.211) we have c^o = 1 
and ci t i = i^ig/2. This yields the following corollary. 

Corollary 4.3. For k = 1, . . . , N , we have 

: S fc (L nr + E)(xi, . . .,x N ) : fCi(x,y) = ^ : t, k (L m + E)(-yt, -y N -i) ■ 

+ ^ : t k _ x (L nr + E) . . . , -y N ^) : J K x (x, y). (4.37) 
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To obtain a counterpart of Corollary 2.4, we take £ = N and note that ( I4.24p implies 

E k (L nT (xi, . . . ,x N ,ifigN/2, . . .,ifigN/2) + E(x u . . .,x N )))C N (x) 

%J f\ S k (l,3,...,2N-l))C N (x). (4.38) 

In this identity we can cancel JC N (x) and the factor (i h ig/2) k , cf. (OljHOB]) . Then it 
says that the matrix 

M = N1 N + Q, Q jk = -5j k coth(z j - z ri ) + (1 - 5,-*) sinh^ - z fc )~\ (4.39) 

has the same symmetric functions as the matrix diag(l, 3, . . . , 2N — 1). Shifting the latter 
matrix and M by iVljy, this leads to the following remarkable result. 

Corollary 4.4. The matrix Q given by (I4.39P has spectrum {— N+l, — iV + 3, . . . , N— 1}. 

We note that, by exploiting the commutation relation 

<p(x, y) : t k (L nr + £) (x) : = : t k \L aT + E + (x) : <^(x, y), (4.40) 

where 

N N-l 



<p(x, y) = exp I — I X m ~ Vn ) ) > ( 4 - 41 ) 

\ \m=l n=l / / 

and the expansion 

A (iiigi\ l (N - k + l\ - w , 
= E( v i ^J ( / J : (2^ + £?) (ar) :, (4.42) 

the exponential factor can be removed from (14.231) . yielding corresponding identities. 

Alternatively, these identities can be deduced from the relativistic kernel identities in 
Theorem 12.11 by using that each PDO : Ti k {L m + E) : can be obtained as a limit of a 
linear combination of the identity operator and the commuting A AOs A±, . . . ,A k that 
result from . . . , A ky+ upon substituting (I3.24[) and 

P = 09- (4.43) 
Indeed, from Section 4.3 in [Ru i94] one can infer 

:t k (L nv + E) :=lini o r fc E(- 1 ) fe+J ( A !Zi)^' A » = L ( 444 ) 

3=0 V 7 
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To illustrate this second method we continue by deducing the pertinent identities for 
= 1. Then we can use Corollary 2.3 to obtain 

E(-l) fc+J (* 1 1) Ai(x!, . . . , xnWx, y) 

3=0 V ' 



3=0 

k 1 fN-n-l\ 

^(-i) fe+ "A n (- 2/1 ,...,- 2/iV _ 1 )E(- 1 ) z N _ k ) s ^y)- ( 4 - 45 ) 

■n— n ;— n \ / 



n=0 «=0 



Hence, for k = N the coefficient of each AAO A n (—y) vanishes. Moreover, for k < N we 
have ('Pascal's triangle') 



N — n\ (N - n - 1\ (N - 1 - n 
N-k)~[ N-k ) ~ \N-l-k 



k = 1,...,N- 1, n = 0,...,Jfe. (4.46) 



Substituting this in (14.45 p . multiplying by /3 fc , and taking the nonrelativistic limit (3 — > 0, 
we arrive at the following result. 

Proposition 4.5. Setting 

N jV-1 

y) = J] II t 2 cosh(/i(x m - y n )/2)]-9/\ (4.47) 

m=l n=l 

we have an eigenfunction identity 

: t N (L m . + E)(x) : ^(x.y) = 0, (4.48) 

and kernel identities 

(: S fc (L nr + S)(a;) : - : S fc (L nr + E){-y) :) /Cx(x, y) = 0, fc = 1, . . . , iV - 1. (4.49) 

For £ > 1, the identities in question become quite unwieldy, and we have not obtained 
them explicitly by either of the above methods. 

The self-duality of the relativistic case is not preserved by the nonrelativistic limit, so 
we proceed to study kernel functions involving the dual variables. To begin with, we need 
to revert to the spectral variables p via (11.28p . To ease the notation we omit the hats, so 
that we wind up with AAOs (cf. ([2.1270 ) 

^ = e n sinh (g (Pro - P „ )/2 ) n «p(=f«mw). (4-50) 

|7|=fcmg7 v " mGi 

Their — > limits yield the dual nonrelativistic AAOs 

A ±fc>nr (p) = e n Pm ~ p : TW n ex p(^ m%j> ( 4 . 5 i) 

I =fc mGY me/ 
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(To be sure, they can also be viewed as the AAOs associated with the relativistic rational 
Calogero-Moser systems.) 

We proceed to obtain the nonrelativistic limit of the dual version of the kernel identities 
in Theorem 2.1. The reparametrized kernel function (I2.32p reads 

nyr G{2tt/h, hfc f3{p m - g n - ifig/2)/fi) 
„„ =1 M G(2n/fi, m P(p m -q n + ing/2)/fi) ' 1 ' ' 

As a preparation for the limit flA.281) . we first use the scale invariance fLA27l) to write 

G(2tt/ fi, hf3; f3z/ fi) = G(l, k; Kz/hfi), k = f3hfi/2n. (4.53) 

Next, we note that we may multiply by a constant and shift the coordinates pj in (14.521) 
without losing the kernel property. A moment's thought then shows that we can invoke 
(1A.28I) to obtain from (I4.52p a nonrelativistic dual kernel function 

\r i \ TT TT ^(p™ ~ Qn)/hfi — g/2h) 

^''l.g rfc-ri/^.W e = o,i,...,N. (4.54) 

Likewise, Theorem 2.1 has the following counterpart. 
Theorem 4.6. For any k £ {1, . . . , N} and r G {+, — } we have 

min(fc^) . . 

A Tkjm (p 1 ,...,p N )jc e (p,q) = ^2 ( .)A-T(k-3),w:(qi,---,qN-i)fci(p,q), ( 4 -55) 

where 

i±m,nr(gi, • • • , q N -t) = 0, 171 > N — £, i ,nr = 1. (4.56) 

From this, the nonrelativistic versions of the three corollaries (12.44R — (12.46H of Theo- 
rem 2.1 will be clear. 

For later use, we add that the nonrelativistic limit of the dual weight function reads 

w (p) = TT r ( i (Pm-p n )/hfi + g/h)r(-i(p m - Pn )/hfi + g/h) 

Km<n<N T (i(p m - p n ) / Hfl)T (-i( Pm - Pn ) / Hfl) 



This follows from ( 1A.28I) in the same way as for the kernel functions. As a check, note 
that the dual nonrelativistic Hamiltonians, 

H±k,nr(P) = ^2 II I ~ Z ) 1 I eX P(^^p,J 

\l\=kmel \ Pm Pn / ma 

nf^^^y 72 , (4.58) 



x 



m&I 
nil 



Pm Pr, 



are then related to the AAOs ^4±fc,nr(p) via 

HlAP) = Wnr(p) 1/2 il,Dr(p)W n r(p)- 1/2 , ±1 = 1, . . . ,N, (4.59) 

as should be the case. 
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4.1.3 The periodic Toda case 

In this section we obtain a kernel identity for the nonrelativistic periodic Toda system via 
its relativistic counterpart ( 12.83(1 . For the same reason as for the elliptic and hyperbolic 
limit transitions, this involves the classical relativistic and nonrelativistic Lax matrices. 
Like in the hyperbolic case, we trade the parameters a + and a_ for 2ir/fi and h/3, resp. 
Choosing also 

V = -ln(^g), g>0, (4.60) 
fi 

it follows from flA.40j) with A = 2 that the relativistic [/-function (I2.7ip satisfies 



\imU(x) = l. (4.61) 

/3-+0 



Moreover, (IA.40I) with A = 1 implies that both kernel functions (12.87(1 and ( 12.88(1 satisfy 



)imS ± (x,y) = 1C(x,y), (4.62) 



where 

N 



K(x, y) = exp ( -| (e^ Xm+1 ~ ym) + e"G"»-*"*)) ) . (4.63) 



h 

m=l 



After the substitutions (I3T24D and flOOl . we obtain from fl2T73D . (I277S1) and (I2T8TT1 



N 



A+ + = J2( 1 + 1 2 e^ Xm+1 - Xm+ihm ) exp(-ihpd Xm ), (4.64) 



m=l 

N 



A{ + = (! + 7 2 e^ Xm - Xm - 1 - i?l/3/2) ) exp(-iH/3d Xm ), (4.65) 



m=l 



N 

2 u,{x m+ x-x m +iKp/2)\ 1/2 



m=l 

x(l + 7 2 e^ m " Xm - 1 - i/l/3/2) ) 1/2 exp(-zW a;m ), (4.66) 

where we have set 

7 = Pm- (4-67) 
Expanding these AAOs in a power series in /3, we get in each of the three cases 

N 



N + /3 [- in ^d x .^+/3 2 H nI + 0(/3 3 ), (4.68) 



where 

k2 N N 
Hnr(x) = — - ^2 C + « 2 E ^^"H a = m (4.69) 



2 

m=l m=l 



is the defining Hamiltonian of the nonrelativistic periodic Toda system. This implies that 
K,{x,y) yields a kernel function for the defining Hamiltonian. Now this can easily be 
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checked directly, but the kernel property for the higher order commuting Hamiltonians is 
harder to show. 

To obtain this more general property, we first need more information on the relation 
between the commuting relativistic AAOs and the commuting nonrelativistic PDOs. To 
this end we begin by recalling a Lax matrix for the classical nonrelativistic periodic Toda 
system. This N x N matrix depends on a spectral parameter w£C* and is given by 

(T \ — X n-l-A _l_ /l 2 A „n( x m-Xm-l) 

\- Lj nT)mn u mnPm "T ^m,n—X ' u "r?i,n+l c 

- {ia) N w8 mN 8 nl - a 2 (za)- JV W - 1 5 ml 5^e^ 1 -^ ) . (4.70) 

Our subsequent considerations involve the symmetric functions of L m , given as the coef- 
ficients Sfc nr in the expansion 

N 

fet{l N + AL nr ) = ^Knr- (4.71) 
fc=0 

In particular, this yields 

N N 

Si, nr = ^p m , S 2 , nr = Yl PmPn-a 2 ^^'^- ( 4 - 72 ) 

m=l l<m<n<N m=l 

(For N = 2 one should add the constant a 2 {w + 1/w) to E 2)nr -) The formula (I4.7ip has 
an unambiguous quantum analog. Indeed, a term in the expansion of the determinant 
that depends on a given p m does not depend on x m , since the only elements of L nr that 
depend on x m occur in the mth column and row. Denoting by L nr the matrix obtained 
from L nr by substituting —ihd Xm for p m , m — 1, . . . , N, the determinant of L nr can now 
be defined by the usual expansion, since no ordering problems arise. Specifically, we have 

N 

det (l N + AL nr ) = Xk ^k,nr, (4.73) 

k=0 

where Sfc,nr is obtained from £fc >nr via the above substitution. Next, we show that the 
PDO S^nr can be obtained as a limit of a certain linear combination of the identity 
operator and commuting AAOs . . . , r £ {+, — }. The latter are obtained from 
the AAOs Af t+ , . . . , A% >+ by substituting (ECMj) and 

2 

rj = - In 7 t ihfi/2, 7 = Pm, (4-74) 
fi 

respectively. Thus we have (cf. f|2773|) . fl2T75|) and fl2T70|) ) 



4 ±} = J2 II ( 1 + 7V ( " m+1 -" m) ) n expi-ihpdvj, k = 1, . . . , N. (4.75) 

\I\=k mel mel 

Lemma 4.7. Fix w E C* , and let 

^ , =(- i ) fc (/^B- i ) fc+ <^:i)^r. w 



k 
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where k — 1, 



iV ; r = +, — , and 



c j 



(1 + {ij) N w) 



N-h 



(1 + (ij^w)"' 1 (1 + (h)^" 1 ) , j = N. 
Then for all k G {1, . . . , iV} and r G {+, — } we have 



(4.77) 



(4.78) 



t R n ~ m h F 1 

- 1 -'mn P u m I -'mm 



Proof. The following reasoning involves a modification of the arguments employed in Sub- 
section 4.3 of [Rui94j to handle the nonrelativistic limit of the elliptic AAOs. We first 
note that in the AAOs D^' all partial derivatives occur to the right of the x-dependent 
coefficients. Since no ordering problems arise in the definition of the differential operators 
Sfc jnr , it suffices to prove the classical version of the statement, obtained by replacing 
-ihV x by p. 

To this end we use the Lax matrix L nr introduced above, as well as a Lax matrix for 
the classical relativistic periodic Toda system. Specifically, let L be the N x N matrix 
given by 

(4.79) 
(4.80) 

(4.81) 
(4.82) 

E ^= ' rn = 2,...,N, (4.83) 

E mn = -(i 1 ) N w, n-m = -2,...,-N+l. (4.84) 



where 
and 



(1 + 7 V (xm+1 -* m) ) i/ *(l + 7 



(/ 2 e n(x m -X m -l)\ 1 / 2 e PP 



Em 
E 



I _|_ ry2 e )l(x 1 -X N ) 

-m = N -2 



1 ? n-m = JV-2,...,l,0, 



By exploiting a limit from a Lax matrix for the elliptic system it can be shown that the 
symmetric functions of L read 



T, k (x,p) = c k S k (x,p), k=l, 
Here, the coefficients c k are given by (14.771) and 



S k (x,p)=J2 II (1 + 7 2 e 



2 u(x m+1 -x m )\ V 2 



r n (i+A 



N. 



fi(x m —x m -i)\ V 2 



(4.85) 



\I\=k rnel 



m-l(£I 



(The details can be found in Subsection 3.2 of [Rui94j .) Moreover, recalling 

7 = Pt*9 = 

one readily checks 

L = l w + /3L nr + 0(/3 2 ), /3^0. 



(4.86) 
(4.87) 
(4.88) 
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Expanding the determinant of the matrix ljv + A/3 1 (L — 1^) in powers of A, we deduce 
the limit 



W',p) = tor 1 ((-!)*(/*) +D-D' +, g:i)^(x,rt 

jfc = l,...,JV. (4.89) 
To relate this limit to A^\x,p), we introduce the diagonal matrix 

(\ _|_ ^,2 e /i(a; m+ i-a' m )\ 1 / 2 

£> = diag(di,...,dtf), rf m = -^ ^ f^. (4.90) 

M + ^/2g/i(a; m -x m _i) j ' 

The principal minor of D with indices {zi, . . . , ik} = I reads 

= n W/l W (4.91) 

so the symmetric functions of the matrices D ±X L are given by 

(D ±X L) (I) = W^CO = ^Af\x lP ). (4.92) 
|/|=* |i|=fc 

Now from f!4.88|) and (14.901) we clearly have 

= ljv + /3L nr + 0(/3 2 ), /3->0. (4.93) 

Therefore, the limit (I4.89P is still valid when we replace Sj by or Aj. Hence the 
classical version of the statement follows, thus completing the proof. □ 

In particular, this lemma implies the commutativity of the N PDOs Si jn r ; • • • > ^N,m 
(a property left open in |Rui94] ) . It is easy to verify that for k — 1 and k = 2 one 
still gets the limit (I4.78P when Ay in (I4.76P is replaced by A^ + , Aj^ + or Hj + with the 
substitutions (13.241) and (14.601) in force. Already for k — 3, however, this is no longer clear, 
since the /3-dependence in the exponentials yields unwieldy extra terms in the expansion, 
cf. e.g. (I4.64p -( )4.66p . (In the elliptic case such extra terms in the AAOs ifi i+ ,if 2 ,+ do 
give nontrivial ^.-dependent deviations from the classical expansion already for k = 2, as 
exhibited in (JUTJ) and (Ojl .) 

We are now prepared for the nonrelativistic counterpart of Theorem 2.7. 

Theorem 4.8. Setting 

D(x) = det (l N + XL nv {x)) , (4.94) 



where L nr (x) is defined by (14.701) with the substitution —ih~d Xm for p m , we have a kernel 
identity 

(D(x)-D(-y))JC(x,y) = 0, (4.95) 

with the kernel function K,(x, y) given by (14.631) . This identity also holds true when /C(x, y) 
is replaced by JC(a(x),y), with o any cyclic permutation. 
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Proof. The limit ( I4.62p is still valid when 77 in the kernel functions S ± (x,y) is replaced 
by ( 14.74ft instead of ( I4.60p . Indeed, this change gives rise to shifts of the arguments of 
the functions Gr and Gl in (I2.87[) ( [2.88j) by ±ihf3/4. This amounts to a shift of z by 
±ia_/4 in the limit formula (1A.40|) . Since the convergence in this formula is uniform on 



compact subsets of the plane, the limit is unchanged. In view of the above lemma and 
(12.831) — (1 2 . 8 4 j) we can now deduce (14.951) . The kernel identities in Theorem 2.7 also hold 
for cyclic transforms of S ,± (x, y), so the same is true for (14.951) . Alternatively, one can 
invoke the cyclic invariance of the PDOs Si iIir , . . . , £jv,nr; which follows from that of the 
AAOs by the lemma. □ 

4.1.4 The nonperiodic Toda case 

As in the periodic case, we first reparametrize a + , a_ and r\ via ( I3.24p and f l4.60j) . Then it 
follows as before (cf. (T4T611 -( l462|) ) that the [/-function ( ET991) has limit 1, whereas 
the nonrelativistic limit of the nonperiodic kernel functions S (x, y)(cf. (12. 1011) — ( 12. 1021) ) 
is given by 



)C(x, y) = exp -| e^ m ~ XN) + ^ (e^ Xm+1 ~ ym) + e^ ym ~ Xm) ) . (4.96) 



m=l 



Proceeding as in the periodic case, we get again (I4.64[) ( |4.67[) with the nonperiodic Toda 
convention ( 11.181) . Hence we deduce the /3-expansion ( 14.681) . now with the defining Hamil- 
tonian 



t2 N 

Hnr(x) = -y E 9 L + a2 E e^~ x - \ a = m (4.97) 

m=l m=l 

of the nonrelativistic nonperiodic Toda system. 

By contrast to the relativistic case, there seems to be no limit transition leading from 
the periodic kernel function given by (14.631) to the nonperiodic one (I4.96p . On the other 
hand, it is again easy to obtain the nonperiodic Hamiltonian (14.971) from the periodic one 
given by (14.691) . It is also straightforward to verify directly the kernel function property 

(H m (x) - H m {y)))C{x, y) = 0. (4.98) 

To obtain the kernel property for the higher order commuting PDOs, we proceed as in 
the periodic case. Thus, we use a nonrelativistic Lax matrix 

{.Lnr)mn &mnPm ^m.n—1 @> &m,n+l& ^ j (4.99) 

and a relativistic one L (14. 79 p . with b m given by 

b x = (1 + 7 2 e^ 2 ^ l) ) 1/2 e /3pi , b N = (1 +1 2 e ^~ x ^) 1/2 (4.100) 
and by (I4.80p for m — 2, . . . , N — 1, and with the matrix E given by 

E mn = l, n-m = N- 1,...,1,0, (4.101) 

E ^= l + ^n-^y m^2,...,N, (4.102) 



57 



E mn = 0, n-m = -2, ...,-N + l. (4.103) 

It is not difficult to check that L m and L are again related by ( 14.881) . and that the 
symmetric functions of L are given by fj4.86j) (with the convention (11.18P in effect). Defin- 
ing the auxiliary A AOs via (14771) and ([4775]) . and then D { k ±] by (14776]) with the 
coefficients Cj set equal to 1, we obtain (14.781) by an easy adaptation of the proof of 
Lemma 4.2. The remark below the proof applies here, too, and now the following analog 
of Theorem 4.3 readily follows. 

Theorem 4.9. Setting 

D{x) = det (l N + XL nv (x)) , (4.104) 



where L m {x) is defined by (I4.99j) with the substitution —ih~d Xm for p m , we have a kernel 
identity 

(D(x)-D(-y))JC(x,y) = 0, (4.105) 
with the kernel function )C(x,y) given by ( 14 . 9 6 1) . 

Finally, we obtain a counterpart of Corollary 2.9, namely, kernel identities that relate 
the PDOs Efc.nr in N variables x = (x\, . . . , Xn) to the PDOs Efc >nr in iV — 1 variables 
V = (l/i, • • -,VN-i) for k < N. 



Corollary 4.10. Setting 



/d(x, y) = exp I -| {e^ Xm+l ' Vm) + e^ ym ~ Xm) ) J , (4.106) 
we have an eigenfunction identity 

Sj V , nr (a;)X:i(x, 3 /) = 0, (4.107) 

and kernel identities 

(S Mr (x) - h^-yfildix, y) = 0, k = I, . . . , N - 1. (4.108) 

Proof. The kernel function /Ci((xi, . . . ,£tv), (yi, . . . , yjv-i)) is obtained from the kernel 
function fC(x,y) given by ( I4.96|) upon substituting y N y N — A and then taking A to 
infinity. Doing so in the Lax matrix L nr (y), the matrix element (Lnr)Ar,JV-i vanishes in 
the limit. Noting d VN annihilates K,\{x,y), the assertions now follow from Theorem 4.4 
upon expansion of the determinants. □ 

4.1.5 The dual nonperiodic Toda case 

To obtain the nonrelativistic limits of the quantities in Subsection 2.5, we proceed in the 
same way as for the dual hyperbolic case. Thus we reparametrize a + ,a_ via (I3.24[) and 
revert to the spectral variables p via (jl.28p . Omitting the hats on p, we wind up with 
A AOs (cf. (127I%2|) and (12TT46]) ) 

= (=*)«"-» e n 2smh(/3(P l- P „) /2 ) n «p(=n>wv.). (4.io9) 

n^7 
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H 



±k, 



En 

\I\=km£l 
n4I 



2sinh(/3(p m -p n )/2) 



1/2 



m&I 



X 



n 

me/ 



related by the weight function 



2sinh(/3(p m -p n )/2) 



1/2 



(4.110) 



W' = ] [ 4sinh(/3(p i - p k )/2) sinh(vr(p j - p k )/hfi), 

i<i<fe<JV 



(4.111) 



cf. (l2~T3gj) - f[2~T39l) . 

Clearly, we need only multiply ( l4TTU9l - fT4TTTU]) by {h^) k{N ~ k) and take p to to get 
the commuting AAOs 



i 



, . r Pm Pit , 
\I\=km£l m£l 

n£I 



H±k,M = e n 



|/|=fcme/ 
n<£I 



Pm Pn 

They are related by f)4.59p . where 



1/2 



JJ exp(=Fi/i/i^ m ) Jl 



mel 



Pm Pn 



1/2 



(4.112) 



(4.113) 



W m {p) = Y[ 2 ((Pi -Pfc)/^) sinh ( 7r (Pi -Pk)/hfi) 

l<j<k<N 



(4.114) 



is the nonrelativistic dual weight function. 

Next, we obtain the nonrelativistic limit of the kernel identities in Subsection 2.5 by 
arguing as in the dual hyperbolic case (cf. the paragraph containing (I4.53P ). using also 
that we may multiply the kernel function by c\ exp(c2 YliPj ~ Qj)) without losing the 
kernel property. 

Theorem 4.11. Letting I G {±l,...,±iV} and o G {±1}, we have the dual kernel 
function identities 

(A^(p) - A-z,„ r (g))/C( P , q y = o, (4.115) 

(H hnT (p) - H^ n M)W D r(p) 1/2 W nv (q) 1/2 ^(P,<lT = 0, (4.116) 

where 



N 



(4.117) 



j,k=i 

Each term of the AAOs occurring in (14.1151) and (14. 1 16[) shifts |/| coordinates by 
±ihfi, so the kernel property is preserved upon multiplication by products of functions 
that are ift/i-antiperiodic. Hence we can derive the kernel property of l/fC(p, q) from that 
of K.(—p, —q) by using the reflection equation of the gamma function, in the form 



T(iz + l/2)T(-iz + 1/2) = n / cosh(Trz). 
We proceed with the following counterpart of Theorem 2.12. 



(4.118) 
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Theorem 4.12. Define kernel functions 

N N-t 

£t(p,q)= H]lT(l-i(p m -q n )/hfi), (4.119) 

m=l n=l 

where £ = 0, 1, . . . , N — 1. For any k 6 {1, . . . , N — £} we have 

A k ,nv(Pi, ■ ■ .,p N )fce(p,q) = ^-fc,nr(?i, • • ■ , q N -i)JC t (p, q) . (4.120) 

Proof. We bypass a laborious derivation from Theorem 2.12 by adapting its proof, as 
follows. First, the analogs of the equations (I2.154p -f l2.155p read 



Ktip, q)' 1 exp(-ihfid Pm )1Ci(p, q) = (ihn) N ~ e TT , (4.121 

AA Pm - q n 



n=l 
N 



/C,(p,g)- 1 exp(^ / i9 gm )/C,(p,g) = (4-122) 

n=l Pn Qm 

Second, we take 5 = + in the identities (12. 153ft and substitute (recall s+(z) = sinh(7T2;/a + )) 

v = ta + p/7i, w = ta + q/-K. (4.123) 

If we now multiply both sides of ( 12. 153ft by t k ( 2N ~ k ~^ and take t — v 0, then we obtain the 
identities 



e n— ^- n — L - 

ic{i,...,N}mei Pn Pm meJ Pm q ' n 

\I\=k np ne{l,...,N-£} 

.^ r ,^^. r „> W'j qm 9n ^"^~ T Pn 9ra 

IC{l,...,N—£} m€l m£l 

\I\=k np ne{l,...,N} 

Finally, using (I4.112p and (j4.121D - fl4.122p . it easily follows that the kernel identities (14.1201) 
amount to fl4.124j> . □ 

Introducing 

K t (p,q) = K: e (-p,-q), (4.125) 

we can mimic the reasoning leading to Corollary 2.13 to obtain the following nonrelativistic 
analog, which concludes this subsection. 

Corollary 4.13. We have eigenfunction identities 

Ae >nT (vi, v N )iCi(v, w) = K, e (v, w), (4.126) 

and kernel identities 

A k ,nr(vi, v N )iCi(v, w) = i_( fc _£), nr (wi , . . . , w N -t)Ki(v, w) , k — £+1, . . . , N. (4.127) 
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4.2 Backlund transformations 



4.2.1 The elliptic case 

To obtain the nonrelativistic versions of the results in Subsection 3.1, we can proceed in 
two distinct ways. First, we can adapt the reasoning based on the expected relation (ll.3ip 
to the nonrelativistic kernel function ^f m (x,y) given by (14.51) . The second way is to take 
j3 to in the relevant formulas in Subsection 3.1. This yields the same results, provided 
we replace p again by i(3g, as we did in Subsection 4.1, cf. (14. IB . Since we want to view 
g as a real coupling constant, it follows from the limits of ( 13 .6p and ( 13 .7p that we wind 
up with purely imaginary momenta pj and qj. Thus, we run into the same problem as 
we had in Subsection 3.1. We improved the situation in the relativistic elliptic setting by 
requiring to be purely imaginary, while keeping p G i(0,a). From this it is clear that 
we can emulate this improvement by replacing p by (3g before taking (3 to 0, keeping g 
real. 

Doing so, we wind up with a generating function 

F m (x,y) = F Wat {x) + F Wni {y) + F K {x,y), (4.128) 

where 

Fw ni {x) = ^ HR{xj-x k + ia/2)R(x 3 -x k -ia/2)), (4.129) 

l<j<k<N 

N 

F K {x,y) =~gJ2 ln (^' - f*))' ( 4 ' 13 °) 

j,k=i 

Obviously, this can also be obtained via f 1 1 . 3 1 j) and ( 14. 5p . provided we replace g by — ig 
in (0~2} and fl43l) . 

In order to show the Backlund property, we use the relation between the elliptic 
relativistic and nonrelativistic Lax matrices from |Rui 94j . The relativistic one is defined 
by 

L it = expfe) jg /(*, - ».) • $^^fy ("3D 

where the function f(z) is given by ( 11.31) . and A G C is a spectral parameter. Its symmetric 
functions are proportional to the Hamiltonians S k (II. ip . Specifically, 

= s(X)- k s(X-p) k - 1 s(X + (k-l)p)S k (x,p), k = l,...,N. (4.132) 

The nonrelativistic Lax matrix is defined by 

Up to a similarity transformation, it coincides with the elliptic Lax matrix introduced 
by Krichever [Kri80j. (At this point a physicist reader might worry about non- matching 
dimensions, inasmuch as p has dimension [momentum], whereas the coupling constant g 
has dimension [position] x [momentum]. But A has dimension [position], and so does s(A). 
Thus the dimensions work out.) 
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Clearly, the Lax matrices L with p = i/3g and L nr are related in the same way as in 
the periodic Toda case, cf. (14.881) . Hence their symmetric functions are related by 

k / AT -1\ 

£ Mr (x, p) = lim p- k ^(-l)** I _ ) E fc (x, p), k = 1, . . . , N. (4.134) 

The same relation holds for the symmetric functions of the matrices L with p = (3g and 
L nr with g replaced by —ig. Therefore, the nonrelativistic Backhand property 

S nr (x,p) = E m (y,q), (4.135) 

follows from its relativistic counterpart (13.81) by using (I4.132p and (14. 134|) . 

Although this reasoning is formally impeccable, it skirts the existence problem already 
discussed in Subsection 3.1. Moreover, the noncompleteness of the flows is now even more 
conspicuous due to the 'wrong sign' of the coupling. Indeed, for the defining Hamiltonian 
associated with L nr (I4.133P we get 



11^ 1 
H m = ^Tr Ll = + 9 2 E P(^-x k )--g 2 N(N-l)p(X). (4.136) 

3=1 l<j<k<N 

Thus, taking g — > —ig leads to a negative coupling, so that the singularities at coinciding 
positions cannot be avoided. (Recall we require that g be real to avoid imaginary momenta 
resulting from the nonrelativistic versions of (13. 6 j) and (13. 7p . cf. (14. 128[) f l4.130[) .) 

4.2.2 The hyperbolic case and its dual 

Proceeding as in the elliptic case, we obtain from (14 . 11 j) — (14 . 1 2 j) once more ( I4.128[) (14. 130[) . 
with R(z) replaced by 2 cosh(pz/2). In the hyperbolic versions of the Lax matrices ( 14.1311) 
and (I4.133P we can take A to infinity. Omitting a similarity factor and reparametrizing 
via (l3T24T) -( l3~25l) and fl4TT]) . we obtain 

/ x . 1 r ( sin 2 (r) \ x l 2 isin(r) 
L r U = exp (Bpi TT 1 + , — ■ — — 4.137 

UJ t£\ smh\p(x J -x l )/2)J smh{ir + »{x j -x k )/2y 1 ; 

L nr (g) jk = 5 jkPj + (1 - M 2sinh(/i( ^_ Xfc)/2r (4-138) 

in accord with |Rui88] (where a parameter z is used instead of ir). Then (I4.88P holds 
true again, so it still holds for the Lax matrices L(—if3pg/2) and L nr (—ig), where we now 
think of (3 being purely imaginary and p, g real. Hence the symmetric functions of the 
latter Lax matrices are again related by (I4.134p . and accordingly the negative coupling 
Backhand property (14.1351) results. As in the relativistic case, its precise interpretation 
within the confines of global analysis/symplectic geometry remains to be determined. 
Turning to the dual case, the relevant dual Lax matrices are [Rui88j 

" / TT ( sin 2 (r) \ isin(r) 
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L Mlt - exp(^) n (i + 1/2 ■ <w .%_ fe) . 

with L nr (g) being the /3 — )■ limit of L(/3fj J g/2). Thus we can expect the nonrelativistic 
Backlund property 

Sk,nr(P,x) = S k ,nr(q,y), (4.141) 

as a limit of the relativistic one. However, it is not straightforward to obtain the nonrel- 
ativistic generating function F nT (p, q) as a limit of F(p, q) (given by (I3.30I) - (I3.32I) with (3 
and \i interchanged), since the renormalizations mentioned below (I4.53P must be taken 
into account. 

Instead, we start directly from the dual kernel function 

^nr(p, q) = W ni (p) 1/2 W ni (q) 1/2 £(P, q), (4-142) 

where /C is defined by (14.541) with i = 0, and W nr by (14.571) . As before, one might expect 
(11.311) to yield the desired generating function F nr (ft, q) as a limit 

limz/iln\I/ nr (p, q). (4.143) 

7lJ,0 

In fact, however, this limit does not exist, as is obvious from the following lemma. 
Lemma 4.14. For ft and q in Gh yP (13.3811 we have classical limits 

~i Ji {Pj-4k)/n+g/2 



limih\n(exp[N 2 (g/h)\n(l/h)]IC(ft,q)) = i V / dwlnw, (4.144) 



ri(pj-Pk)/n+g 

liminin(exp[-JV(iV- l){g/h)\n(l/h))W nT (p)) =i V / dwlnw, (4.145) 

m j<k Ji(p-j-p k )/»-g 

where the integration paths stay away from the cut (— oo, 0]. 
Proof. We recall Stokes' formula, in the form 

lim -j-ln (exp[(rf-M)AlnAi r ^^ ) = f dwlnw, u, d £ (-oo, 0]. (4.146) 
A^ooA ^ T(Ad)J J d 

Inspecting the definitions (14.541) (with i = 0) and (I4.57p . the limits readily follow from 
this. □ 

Instead of using ( 14. 143ft . it is now clear that we need to define the generating function 

by 

F m (ft, q) = \imih\n(exp[N(g/H) ln(l/ft)]^ nr (p, q)). (4.147) 

Hence q(x,ft) is to be determined from the equations 

dF ni 



Xj 



J 3 

N 



-Yln( Pi-P'-W) + 1 -yin( Pj^h+Ml) , (4.148) 
2 V f£ \Pj -Pk + ^gJ ^ ~ \Pj - Qk - W 2 / 
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and then y(x,p) is given by 
dF nv 



= -yin(j j ~ ft + ifl9 ) +lyin( Pk -^ + ^f] . (4.149) 

Clearly, when we retain the physical choice fi, g > 0, then we get a contradiction from 
assuming that for given (x,p) G Cl hyp the implicit equations (j4.148j) yield a solution 
g G Ghyp. (Indeed, it would follow that x is purely imaginary.) We can only avoid this 
snag by keeping \i positive, while requiring that g be purely imaginary. 

Accepting this and assuming (possibly complex) solutions, we deduce as before (cf. Sub- 
section 3.1) that the Backhand property is equivalent to the functional identities 

y> 1 r Pm-pn- I r Pm - q>n + iyg/2 

/C{WV>mA Pm ~ Pn M Pm-qn-iM/2 

\I\=k n n=l,...,N 

_ TT Qm ~ q n + ifig tt Pn - q m + ifig/2 ^ 

|J| = fc n^-f n=l,...,A r 

These are easily deduced from (13.33!) . an d so f )4.14ip follows. 
4.2.3 The periodic Toda case 

Comparing the expected asymptotic relation (jl.3ip to the nonrelativistic kernel function 
K,{x, y) given by ( 14.63[) . it becomes clear that in this case no limit is needed. The resulting 
generating function, 

N 

-igY^ (e^ Xm+1 - ym) + e» {ym - Xm) ) , (4.151) 

m=l 

would yield purely imaginary momenta, but like in Subsection 3.3 this disease can be 
cured by an analytic continuation 

x m — > x m — iir/2/j,, m—l,...,N. (4.152) 

Then we obtain the generating function 

N 

Fnr(x, y)=gY^ {-e^ Xm+1 - ym) + e ^ Vm ~ Xm) ) , (4.153) 



m=l 



which gives rise to 

p m = M {e^ Vm ' Xm) + e^^-w™-!) , (4.154) 
q m = fig ( e ^ Vm - Xm) + e^ Xm+1 ~ ym ) . (4.155) 

These equations can be regarded as the nonrelativistic limit of the equations (I3.5ip - p.52p . 
Indeed, replacing 7 by (3 fig and shifting x m — > x m — \n((3fig)/ fi in the latter equations (in 
accord with <K0T\\ and (J32ZD)> & is obvious that their (3 -> limit yields fl4.154p - fl4.155D . 
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Next, we recall that we already detailed a Lax matrix L nr (I4.70P for the nonrelativistic 
periodic Toda system, and obtained its symmetric functions £fc inr as limits of appropriate 
linear combinations of the relativistic Hamiltonians S k , cf . ( 14. 89ft . Therefore the Backlund 
property 

£ fcinr (x,p) = E k;ia (y,q), k = l,...,N, (4.156) 
follows from its relativistic counterpart. For the defining Hamiltonian 

N N 

H m {x lP ) = -Y,P 2 rn + * 2 Y. e " {Xm+1 ~ Xm) ' a = M, (4-157) 

771=1 771=1 

it is of course easily checked directly from (I4.154l) - (l4.155p . 



4.2.4 The nonperiodic Toda case and its dual 

The nonperiodic Toda counterparts of the formulas (14 . 1 5 1 j) (14 . 1 5 71) will be obvious by 
now: we need only insist on the convention (11.181) for x and y. Also, the Backlund 
property follows by using the Lax matrices given by (I4.99j) — (14. 103h . 

A study of the dual system is less straightforward. Of course, from ( 14.1141) we get the 
same result 

limiMnW„ r (p) = — V (p m -p n ), (4.158) 
HO U — ' 

^ \<m<n<N 

as in the relativistic case, cf. (I3.63p . For the kernel functions lC(p,q) cr given by (14. 1 lTj) . 
however, the relevant limit vanishes for a = — 1 and does not exist for a = 1. On the other 
hand, in this case we have even more freedom to modify kernel functions than indicated 
in the paragraph containing ( ll.29j) . since we can also multiply by products of arbitrary 
ifo/z-antiperiodic functions. Exploiting this, one can probably obtain the same generating 
function as we now shall arrive at by starting from the modified relativistic generating 
function F(p, q) (ETTOl . 

The crux is that when we add the function 

N " 

--\n(/3)J2(Pm-q m ), (4.159) 

771=1 

to F, discard a constant, and then take (3 to 0, we get the limit function 
F m (p,q) = y J2 {(N-2m+2)p m +(N-2m)q m j + - J2 / dwln(-iw). (4.160) 

^ 777=1 ^ 771,77=1 " 

As in Subsection 3.4, the corresponding equations 

N 



OF, 



~^{N - 2m + 2) - - J^ln ( - i(p m - q n )), (4.161) 



dp m 2/i II 
dF nr in 



77=1 

N 



(JV _ 2m) - - J2 ln ( - l (Pn ~ Qm)) , (4-162) 



dq m 2/i n n t 
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might well yield a solution q e G, y G K^, with the interlacing property (I3.72p . In any 
case, the Backhand property 

H k ,m{x,p) = H kjDI (y,q), k = l,...,N, (4.163) 

for the dual Hamiltonians 

H k!m (x,p)= Yl II I - - I II ex P(^)' k = l,...,N, (4.164) 

icz{x,...,N}mei * Pm Pn * lei 

\I\=k nt 1 

can now be shown in the same way as in Subsection 3.4. 



A Elliptic and hyperbolic gamma functions 

The elliptic and hyperbolic gamma functions were introduced and studied in [Rui97j as 
so-called minimal solutions of certain first order analytic difference equations. (The hy- 
perbolic gamma function has various differently-named cousins, as detailed in Appendix A 
of |Rui05j .) In this appendix we review features of these gamma functions that are relevant 
for the present paper. 

The following material concerning the elliptic gamma function G(r,a + ,a_; z) can all 
be found in Subsection III B of [Rui97j . To begin with, the elliptic gamma function can 
be defined by the product representation 

n°° 1 — exp ( — (2m + l)ra+ — (2n + l)ra_ — 2irz) 
—, — f- (A- 1 ) 
m n _ 1 — ex P ( — (2 m + l) ra + — (2w + l)ra_ + 2irz) 

Here and below, we require that the parameters satisfy 

r, a+, a- > 0. (A.2) 

It is obvious from ( 1A.ll) that the elliptic gamma function is meromorphic in z, with poles 
and zeros that can be read off. In particular, for z in the strip 

|S(s)|<a, a=(a + + a_)/2, (A.3) 

no poles and zeros occur, so that we have 

G(z)=exp(ig(z)), |3(z)| < a, (A.4) 

with the function g(z) being analytic in the strip. (We often suppress parameters when 
no ambiguity arises.) In fact, it is explicitly given by 

/ \ sm(2nrz) Nl , . 

9 {r,a + ,a_;z = > - — — — — -, \$S(z) \< a. A.5 

2nsmh{nra + ) smh(nra_J 

Both from this series representation and from (I A . 1 j) . the following properties are clear: 

G(—z) = 1/G(z), (reflection equation), (A. 6) 
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G(a_, a + ; z) = G(a + , a_; z), (modular invariance), (A. 7) 

G^^r, \a + , \a^; Xz) = G(r,a + ,a_; z), Ag(0,oo), (scale invariance). (A. 8) 

The elliptic gamma function arises as a minimal solution of analytic difference equa- 
tions that involve a right-hand side function defined by 

oo 

R(r, a; z) = J]jl - exp{2irz - (2k - l)ar)} [1 - exp(-2irz - (2k - l)ar)]. (A.9) 

k=l 

Specifically, setting 

R 5 (z) = R(r,a 5 ;z), 5 = +,-, (A.10) 

it satisfies 

G(z + ia s /2) 

GH.,ffl ° 6 -' a ' = +■- (A - n) 

It is clear from (IA.9j) and ( lA.lOp that the functions R± are entire, even and 7r/r-periodic, 
and satisfy 

R,(; + m /2) f= (A12) 

R 5 (z - ia s /2) 

For the classical and nonrelativistic limits in the main text we need to invoke two 
related zero step size limits of the elliptic gamma function. First, for z and w staying 
away from cuts given by 

±i[a/2,oo) + kn/r, k G Z, (A.13) 

we have z 

lim a_^(r, a, a_; z) = — I dwln R(r, a;w), (A. 14) 

where the logarithm takes real values for w real; moreover, this limit is uniform on compact 
subsets of the cut plane. In particular, when we have an upper limit z with ^(z)! < a/2, 
we can choose a path along which < a/2 and use the representation 

lnfl(r,a;z) = -f; cos ( 2 ^) |^)| < a/2, (A.15) 

n=l v ' 

which follows from ( lA.llj) and ( 1A.5j) . The second limit reads 

G(r, a, a_; z + iuaJ) 

lim— — r = exp((M — a) hxRir, a; z)), u, a E K, (A. 16) 

a-4-o G(r, a, a_; z + idaJ) u ' v " v ' 

uniformly on compact subsets of the cut plane. Note that for u — d integer this limit 
readily follows from (lA.llj) . 

In the main text we make extensive use of the functions 

s s (z) = s(r, a s ;z), 5 = +,-, (A. 17) 

defined by 

s(r, a; z) = exp(—r]rz 2 /ii)a(z; 7i/2r, ia/2), (A. 18) 
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where a denotes the Weierstrass sigma function. Hence these functions are entire, odd, 
7r/r-antiperiodic and satisfy 

s s (z - ia s /2) 
They are related to R± via the formula 

1 00 1 

ss{z) = 77- TT ~, 3 " exp(irz)Rg(z + ia s /2). (A.20) 

2ir t=i (1 - exp(-2ka s r)) 

Also, the well-known product representation for the Weierstrass sigma function entails 

f\ ^ ( 2 / \ • 1/ / X TT (j~ ~ 2?r/ / r )) ( Z ~» ~g) /A 01 N 

s*(z) = — e 5 (-r2; /vr) smh(7rz/a 5 ) I I -^—2 (A.21) 

where we have introduced the notation 

e s (z) = exp(nz/a s ), 6 = +,—. (A. 22) 

Clearly, this implies 

lim sg(z) = — smh.(irz/as) } (A. 23) 

the limit being uniform on compact subsets of C. 

We proceed to discuss the hyperbolic gamma function G(a + , a_; z), cf. Subsection III A 
of [Rui97j . It can be defined as the unique minimal solution of the analytic difference 
equations 

w^W>^ 2mshM - ^ + '"' (A - 24) 

satisfying G(0) = 1. It arises from the elliptic gamma function by the following limit: 

( k 2 z \ 

limG(r, a+, a_; z) exp I = G(a + , a_; z). (A.25) 

rio \mra + a_ J 

It is meromorphic in z, and for z in the strip flA.3[) it has neither poles nor zeros. Thus it 
can be written as in (1A.4j) . with g(z) analytic in the strip. Explicitly, g(z) has the integral 
representation 

g(a +J a_;z) = f°° *L f ^y* . --?—), \*(z)\ < a. (A.26) 

J y \2 smh(a + y) smh(a-y) a + a_y J 

From this it is clear that the hyperbolic gamma function also satisfies the reflection 
equation (1A.6I) and has the modular invariance property (1A.7jl . whereas the counterpart 
of (m> reads 

G(Xa + , Aa_; Xz) = G(a + , a_; z), A £ (0, 00), (scale invariance). (A. 27) 

In the main text we need several zero step size limits of the hyperbolic gamma function. 
The first one yields the relation to the Euler gamma function: 

limG(l,K;K2 + V2)exp(^ln(2vrfi;)-ln(27r)/2) = l/r(iz + 1/2). (A.28) 
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The second and third one are needed for the classical and nonrelativistic limits, resp.: For 
z and w staying away from cuts given by ±z[a/2, oo), we have 

lim a_g(a, a_; z) = — \ dw ln(2 cosh(7ru;/a)), (A.29) 

lim ffi ' °~' Z + lU 7 a ~] = exp((u - d) ln(2 cosh(7rz/a))), u,d6l, (A.30) 

uniformly on compact subsets of the cut plane. 

In the relativistic Toda setting it is expedient to switch to two slightly different hy- 
perbolic gamma functions given by 

> (A.31) 

2a + a_ / 

(%TIZ^ \ 
-i>X ~ n ' (A.32) 

2a + a_ / 

where 

* = k (- + -) ■ (A33) 

24 \a_ a+ / 

These functions are the unique minimal solutions of the analytic difference equations 

G R (z + ia- S /2) 



G R (z - ia-s/2) 

G L (z + ia_ 5 /2) 
G L (z - ia- S /2) 

with asymptotic behavior 



l + e s (-2z), (A.34) 
l + e s (2z), (A.35) 



Gr(z) — l + O (exp(—r\U(z)\)) , U(z) ±oo, (A.36) 

Gfl(z) = exp(±i (2x + vr2 2 /a+a-) )(1 + C(exp(-r |3?(^)|))), -> Too, (A.37) 

where the decay rate r can be any positive number satisfying 

r < 27rmin(a + , a_)/a + a_. (A. 38) 

Furthermore, they are related by 

G*(*)G L (-z) = 1. (A.39) 

The properties of the functions Gr and Gl just stated are easy to infer from the corre- 
sponding properties of the hyperbolic gamma function. See also Appendix A in [Rui05j . 
where functions Sr and Sl were introduced that differ from Gr and Gl by the shift 
z — > z — ia. 

Finally, we have occasion to use the limits 

limG K (a + , ^±A?±l„J-WfP (±^(^)) • * = J. (A.40) 
a„^o l V 2ir o./ \1, A > 1, v ; 

which hold uniformly on compact subsets of C. These limits are proved in Appendix B 
of [Ruillj . 
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B A family of Hilbert-Schmidt operators 

In this appendix we reconsider the periodic Toda kernel functions given by (12.87D and 
(12.881) . More precisely, shifting y n by £ in S + (x,y), we focus on the resulting function 

noting that basically the same results apply to S~(x,y). Letting £ vary over the strip 

|3(0I < a/2, (B.2) 

we shall associate a family of Hilbert-Schmidt operators to the kernel functions (IB. ID . To 
this end we need some preliminaries. 
We denote by E the Euclidean space 

E = {x e R N | x x + • •• + x N = 0}, (B.3) 

with inner product given by the restriction of the standard inner product on R . Fur- 
thermore, we denote the Lebesgue measure on E by A# and the corresponding Hilbert 
space by L 2 (E). Next, we introduce a change of coordinates 

s = — (ari H Kxat), r n = x n -x n+1 , n = 1, . . . , JV - 1, (B.4) 



with inverse given by 

A — J V-l 



a; m = s - ^ ^nr n + ^ r n , m = 1, . . . , iV. (B.5) 



71=1 



The associated Jacobian determinant equals 1, as is readily verified. Viewing r%, . . . , tat_i 
as coordinates on E by taking s = in flB.51) . we deduce 

d\ E = N- 1/2 d ri ---dr N ^. (B.6) 

To explain the relevance of these coordinates for the Toda A AOs A.fg(x) (cf. ( 12. 73ft - 
(I2.74p ). we point out that they factorise as a product of a center-of-mass operator and a 
reduced operator: 

At & = AlrK^ At/ m = exp(-da„ 5 d s /N). (B.7) 

Here, the reduced A AOs A^g depend only on the variables r 1; . . . , r N _i, so they commute 
with the center-of-mass operators. It is easy to verify that the kernel identity (I2.83j) 
remains valid if we substitute for A± IS either A±f s or *A^ C ™. Hence, viewing L 2 ($L N ) as 
a tensor product 

L 2 (R N )~L 2 (R)(g)L 2 (E), (B.8) 

the findings of this appendix can be used to study Hilbert space aspects of the AAOs, 
but this is beyond our present scope. 

We are now prepared to state the main result of this appendix. 
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Theorem B.l. The operator on L 2 (E) defined by 

(4/)(x)= / Ss(x,y)f(y)d\ E (y), f G L 2 (E), (B.9) 

J E 

with S^(x,y) given by ( IB. II) . is a Hilbert- Schmidt operator, provided I; satisfies ( IB. 21) . 



To prove this theorem we need to show that the ^-restriction entails 

h= I \Ss(x, y)\ 2 d\ E (x)d\ E (y) <oc. (B.10) 

JExE 

As it stands, this integral is hard to estimate. Using the above coordinates r on E(x) and 
the same coordinates q on E(y) we have from (IB.6j) 



7 € = N- 1 [ \S ( (x( r ),y(q))\ 2 drdq, (B.ll) 

JExE 

but at face value this seems no improvement. 

We shall therefore introduce new coordinates on E x E, in which the function S^(x, y) 
takes a particularly simple form. First, we set 

U 2n -1 = X n - y n , U 2n = yn-X n+1 , 71 = 1, . . . , N - 1. (B.12) 

Now on E x E we have xi + • ■ ■ + xn = yi + ■ • • + yN = 0, so that 

UN - Xl = -U2 - U4 U 2 N-2, (B.13) 

y N - %n = u 1 + u 3 J < hM 2 jv- 3 - (B- 14 ) 

For (x, y) G E x E, we thus have 

G R (-u 2 —u 4 u 2 n-2 - - r}/2 + £) 



GlK + u 3 + ■ ■ ■ + u 2N „ 3 + ia/2 + 77/2 + 

x n GR ( U2n ~ ia / 2 ~v/ 2 +o 



(B.15) 



A G^-u^ + ia/2 + T]/2 + ' 



With a slight abuse of notation, we shall write S^(u) for the right-hand side of this 
expression. 

To check that ui, . . . ,U2N-2 yield well-defined coordinates on E x E, we first note 
that (lB~T2l entails 

r n = x n - x n+1 = u 2n -i + u 2n , n = l,...,iV-l, (B.16) 

q n = y n ~ Vn+i = u 2n + u 2n +i, n = 1, . . . ,N - 2, (B.17) 

?iv-i = yjv-i - Vn = u 2N ^ 2 — (u 1 + u 3 -\ h u 2N _ 3 ). (B.18) 

From this we readily deduce 

d(r 1 ,q 1 , . . . ,rjv-i,gjv-i) 



(9(«i, . . . ,U 2 N~2) 
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A/", (B.19) 



so the linear transformation given by ( IB.16j) - (lB.18l) maps R 2N ~ 2 onto R 2N ~ 2 . Further- 



more, from (IB. Ill) we get 

Jf= / \S^u)\ 2 du. (B.20) 

Therefore, the inequality (IB. 10|) . and hence the validity of Theorem IB.lt are clear from 
the following lemma. 

Lemma B.2. Assuming £ satisfies ( IB. 21) . there exist constants B, b > sitc/i i/iat 

\3 ( {u)\ < Bexp(-b\\u\\), VnGl 2Ar " 2 . (B.21) 
In order to prove this lemma we shall make use of the following observation. 

Lemma B.3. Let a = (ai, . . . ,a n ) G R n be such that \ \a\ \ = 1. Le£ P and iV fre £/ie swrn 
of the positive and negative aj, respectively. Then, at least one of the numbers P + N and 
—N is greater than or equal to l/(2y/n). 

Proof. If N < — l/(2y / n), then the statement is obviously true. Next, assume N is greater 
than —l/(2s/n). Clearly, \ctj\ > l/\/n for at least one j. By the assumption on N, all 
such oij must be positive. It follows that P > 1/ ^Jn, so that 

P + N>-^= 1 — = -^=. (B.22) 

^ 2^E 2^i v ; 

□ 

Proof of Lemma \B. M We observe that S^(u) is a meromorphic function of u whose poles 
are located at 

U2n = v/2 ~~ £ ~ — i{ka + + Za_), u 2n -i = f]/2 + £ — ia/2 — i(ka + + la-), (B.23) 

with n = 1, . . . , N — 1, and at 

u 2 + Ui + ■ ■ ■ + u 2N -2 = -??/2 + £ + ia/2 + i{ka+ + Zo_), (B.24) 
Mi + m 3 H h M2Af-3 = —»7/2 - £ + «a/ 2 + i(£;a + + Za_), (B.25) 

where fe, Z G N. Since we have —a/2 < < a/2, it follows that S^(u) is real-analytic 
for all u G R 2Ar_2 , and hence bounded on any compact subset of M 27V ~ 2 . 

We proceed to study the asymptotic behavior of S^{u). To this end we assume ti^O, 
and let (3 = (/3i, . . . , P2N-2) be the corresponding unit vector m/||m||. We observe that 
the numerator and denominator in (lB.15j) depend only on (3 n with n even and odd, 
respectively. It is therefore convenient to introduce 

0° = (pi,P3,...,Pw-a), (B.26) 
(3 e = ((3 2 ,f3,,...,(3 2N - 2 ). (B.27) 

It follows from Lemma IB.3I and (1A.36|) (IA.37h that the numerator of (IB.15P is bounded 
above by 

Cexp [ — — \\u\\ ) (B.28) 

a + a„VjV — 1 
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for some constant C > 0. Similarly, the denominator is bounded below by 



a i a. 



for some constant D > 0. This clearly implies the statement. In fact, these bounds imply 
that we can choose 

b = ?=== min(a/2 - 9(0, + a/2). (B.30) 

a + a_VA — 1 

□ 
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